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Abstract 



We investigate the structure of neutron stars shortly after they are born, when the entropy 
per baryon is of order 1 or 2 and neutrinos are trapped on dynamical timescales. We find 
that the structure depends more sensitively on the composition of the star than on its 
entropy, and that the number of trapped neutrinos play an important role in determining 
the composition. Since the structure is chiefly determined by the pressure of the strongly 
interacting constituents and the nature of the strong interactions is poorly understood at 
high density, we consider several models of dense matter, including matter with strangeness- 
rich hyperons, a kaon condensate and quark matter. 

In all cases, the thermal effects for an entropy per baryon of order 2 or less are small 
when considering the maximum neutron star mass. Neutrino trapping, however, signifi- 
cantly changes the maximum mass due to the abundance of electrons. When matter is 
allowed to contain only nucleons and leptons, trapping decreases the maximum mass by an 
amount comparable to, but somewhat larger than, the increase due to finite entropy. When 
matter is allowed to contain strongly interacting negatively charged particles, in the form 
of strange baryons, a kaon condensate, or quarks, trapping instead results in an increase 
in the maximum mass, which adds to the effects of finite entropy. A net increase of order 
O.2M occurs. 

The presence of negatively-charged particles has two major implications for the neutrino 
signature of gravitational collapse supernovae. First, the value of the maximum mass will 
decrease during the early evolution of a neutron star as it loses trapped neutrinos, so that if a 
black hole forms, it either does so immediately after the bounce (accretion being completed 
in a second or two) or it is delayed for a neutrino diffusion timescale of ~ 10 s. The latter 
case is most likely if the maximum mass of the hot star with trapped neutrinos is near 
1.5M . In the absence of negatively-charged hadrons, black hole formation would be due 
to accretion and therefore is likely to occur only immediately after bounce. Second, the 
appearance of hadronic negative charges results in a general softening of the equation of 
state that may be observable in the neutrino luminosities and average energies. Further, 
these additional negative charges decrease the electron fraction and may be observed in the 
relative excess of electron neutrinos compared to other neutrinos. 
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1 A neutron star is born 



A protoneutron star is formed in the aftermath of the gravitational collapse of the core of a 
massive star. Its evolution proceeds through several distinct stages [El, 0], which may have 
various outcomes, as shown schematically in Fig. 1. 

1. Immediately following core bounce and the passage of a shock through the outer 
protoneutron star's mantle, the star contains an unshocked, low entropy core of 0.5 — 
0.7 Mq in which neutrinos are trapped ||. This is surrounded by a low density, high 
entropy mantle that is both accreting matter falling through the shock and rapidly 
losing energy due to beta decays and neutrino emission. The shock is momentarily 
stationary prior to an eventual explosion. 

2. On a timescale of about 0.5 s, accretion becomes much less important as the supernova 
explodes and the shock lifts off the stellar envelope. Extensive neutrino losses and 
deleptonization of the mantle lead to the loss of lepton pressure and to collapse of the 
mantle on the same timescale. Neutrino diffusion times from the core are too long 
to significantly alter the core during this stage. If enough accretion occurs, and the 
initial core were large enough, the mass of the hot, lepton-rich matter could exceed 
the maximum mass which is stable, in which case the star would collapse to form a 
black hole. In this event, neutrino emission would effectively cease, since the event 
horizon is believed to form outside the neutrino photosphere |5|. 

3. This stage is dominated by neutrino diffusion causing deleptonization and heating of 
the core. Neutrino-nucleon absorption reactions set the diffusion timescale to about 
10-15 s. The maximum entropy per baryon reached in the core is about 2 (in units of 
Boltzmann's constant). When the core deleptonizes, the threshold for the appearance 
of strangeness, in the form of hyperons, a Bose kaon condensate, or quarks, will be 
reduced || [?], ^, ^J. If one (or more) of these additional components is present, the 
equation of state will soften, leading to a decrease in the maximum mass. There is, 
therefore, the possibility that a black hole could form at this later time. 

4. Following deleptonization, the star has a high entropy, so that thermally produced 
neutrino pairs of all flavors are abundant, and thermal diffusion and cooling of the 
hot neutron star takes place. Because the entropy is higher at the beginning of 
cooling than it is at the beginning of deleptonization, the neutrino mean free paths 
are smaller and the timescales longer. In approximately 50 s, as the average neutrino 
energy decreases, the star becomes essentially transparent to neutrinos and the core 
achieves a cold, catalyzed configuration. The loss of thermal energy leads to a small 
increase in the threshold density for the appearance of strange matter, so that, in the 
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absence of further accretion, it is unlikely that a black hole could form during this or 
later phases. 



5. Following the onset of neutrino transparency, the core continues to cool by neutrino 
emission, but the star's crust cools less because of its lower neutrino emissivity. The 
crust acts as an insulating blanket which prevents the star from coming to complete 
thermal equilibrium and keeps the surface relatively warm (T ~ 3 x 10 6 K) for up to 
100 years. This timescale is primarily sensitive to the neutron star's radius and the 
thermal conductivity of the mantle p0| . 



6. Ultimately, the star achieves thermal equilibrium when the energy stored in the crust 
is depleted. The temperature to which the surface now cools is determined by the rate 
of neutrino emission in the star's core. If this rate is large, the surface temperature 
will become relatively small, and the photon luminosity may become virtually unde- 
tectable from the Earth. This will be the case if the direct Urea (beta decay) process 
can occur, which happens if the nuclear symmetry energy is large or if hyperons, a 
Bose condensate, or quarks are present. Somewhat higher surface temperatures occur 
if superfluidity in the core cuts off the direct Urea rate below the superfluid's critical 
temperature. A relatively high surface temperature will persist if the Urea process 
can only occur indirectly with the participation of a spectator nucleon - the so-called 
standard cooling scenario. 

Neutrino observations from a galactic supernova will shed much light on the first four 
of the above stages. Observations of X-rays or 7— rays from very young neutron stars are 
crucial for the last stages. The duration of each of these stages is essentially determined 
by neutrino diffusion timescales, and thus depends both upon the microphysics and the 
macrophysical structure of neutron stars. Roughly, the diffusion timescale is proportional 
to i? 2 (cA) _1 , where R is the star's radius and A is the effective neutrino mean free path. 
Thus, important constraints upon the properties of dense matter can be achieved by looking 
at this relation as it applies to each stage. Generally, the structure (i.e., mass, radius, 
etc.) of both hot and cold, and both neutrino-rich and neutrino-poor, stars is fixed by the 
equation of state (EOS) and the composition. Both are also crucial to knowledge of the 
neutrino mean free paths. 

The behavior of the maximum mass as a function of temperature and neutrino trapping 
is of practical importance If the maximum mass of a cold, neutrino- free neutron star is 



near the largest measured masses of neutron stars [12], namely 1.442±0.001 Mq in the case 
of PSR 1913+16, the question arises as to whether or not a larger mass can be stabilized 
during the preceding evolution of a star from a hot, lepton-rich state. If it cannot, then any 
transition from a neutron star to a black hole in the stellar collapse process should occur 
extremely early, either immediately or during the first few tenths of a second, when the 
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proto-neutron star is rapidly accreting mass from unejected matter behind the shock. If a 
larger mass can be stabilized, however, the transition from a neutron star to a black hole, if 
it occurs at all, could happen later, on the neutrino diffusion or thermal timescale, namely 
~ 10 seconds. Burrows || has demonstrated that the appearance of a black hole should 
be accompanied by a dramatic cessation of the neutrino signal (since the event horizon 
invariably forms outside the neutrinosphere). 

A newly-formed neutron star should accrete its final baryon mass within a second or 
two of its birth, so that the neutrinos will not have had time to diffuse from the stellar 
core. Thus, both the maximum mass for a hot, neutrino-trapped star and a cold, catalyzed 
star must be greater than the largest measured neutron star mass. Moreover, due to the 
binding energy released because of neutrino emission and cooling, the maximum mass for a 
hot, neutrino-trapped star must be at least 0.1 — 0.2 Mq larger than this limit. This feature 
could provide a more severe constraint upon the equation of state than limits based solely 
upon cold, catalyzed matter. This is especially true if the neutrino-trapped maximum mass 
is less than the cold, catalyzed maximum mass for a given equation of state. 

2 Scope of this work 

The evolutionary scenario presented above is based on dynamical calculations carried 
out with a schematic equation of state, since little work has been carried out on the struc- 
ture of protoneutron stars shortly after their birth, although, of course, there have been 
many investigations of cold neutron stars. The purpose of this work is to investigate the 
composition and structure of these newly born stars. It would then be of interest to study 
the implications for the dynamical evolution of neutron stars, but this we defer to the 
future. 

There are two new effects to be considered for a newborn star. Firstly, thermal ef- 
fects which result in an approximately uniform entropy /baryon of 1-2 across the star m. 
Secondly, the fact that neutrinos are trapped in the star, which means that the neutrino 
chemical potential is non-zero and this alters the chemical equilibrium, which leads to com- 
positional changes. Both effects may result in observable consequences in the neutrino 
signature from a supernova and may also play an important role in determining whether 
or not a given supernova ultimately produces a cold neutron star or a black hole. 

Since the composition of a neutron star chiefly depends upon the nature of the strong 
interactions, which are not well understood in dense matter, we shall investigate many of 
the possible models. After a brief discussion of the equilibrium conditions in section 3, 
we begin section 4 by discussing non-relativistic potential models and their predictions for 
protoneutron stars. We then turn to relativistic models, which may be more appropriate, 
since the central densities involved are high. At first, we allow only nucleons to be present 
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in addition to the leptons. Neutrino trapping increases the electron chemical potential 
and, therefore, the lepton and (due to charge neutrality) the proton abundances. With 
more protons, the equation of state is softer and the maximum neutron star mass is lower. 
However, it has been recently realized that if additional negatively charged particles, such as 
kaons, hyperons or quarks, are present, this situation can be qualitatively changed || |7], ||, g]. 
This is due to the change in the chemical potentials, which delays the appearance of these 
strongly interacting particles that lower the pressure. Their softening effect is therefore 
less in evidence, and a larger maximum mass is obtained for the young star. This means 
that the star could become unstable when the initial neutrino population has departed, as 
we mentioned in section 1. In the remainder of this section, we therefore investigate in 
some detail the effect of hyperons, kaon condensates, and quarks upon the structure of the 
protoneutron star. 

The evolution of the protoneutron star is determined by the time scales involved. While 
accurate results require detailed numerical work that incorporates neutrino transport, it is 
nevertheless possible to use semi-analytical techniques to gain an understanding of the times 
involved and the general energetics. This is discussed in section 5. 

In section 6, we discuss the implications of this work for delayed black hole formation 
and for neutrino signals from supernovae, in general, and SN1987A in particular. 



3 Equilibrium conditions 

For stars in which the strongly interacting particles are only baryons, the composition 
is determined by the requirements of charge neutrality and equilibrium under the weak 
processes 

B 2 +£ + Vj,; B 2 +£^Bt+u e: (1) 

where B\ and B 2 are baryons, and I is a lepton, either an electron or a muon. Under 
conditions when the neutrinos have left the system, these two requirements imply that the 
relations 

EH^+nW) = £(n« + nH) (2) 

i i 

Hi = hfin - qifie , (3) 

are satisfied. Above, n denotes the number density and the superscripts (±) on n signify 
positive or negative charge. The symbol fii refers to the chemical potential of baryon i, bi 
is its baryon number and qi is its charge. The chemical potential of the neutron is denoted 
by fi n - 
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Under conditions when the neutrinos are trapped in the system, the beta equilibrium 
condition Eq. (0) is altered to 



(4) 



where \i Vl is the chemical potential of the neutrino U£. Because of trapping, the numbers of 
leptons per baryon of each flavor of neutrino, i = e and [/,, 



Y U = Yi + Y t 



VI 1 



(5) 



are conserved on dynamical time scales. Gravitational collapse calculations of the white- 
dwarf core of massive stars indicate that at the onset of trapping, the electron lepton 
number Yi e = Y e + Y Ve ~ 0.4, the precise value depending on the efficiency of electron 
capture reactions during the initial collapse stage. Also, because no muons are present 
when neutrinos become trapped, the constraint Yl^ = Y/j, + Y u = can be imposed. We 
fix Yl£ at these values in our calculations for neutrino trapped matter. 

For completeness, we give here the partition function for the leptons. Since their 
interactions give negligible contributions [13], it is sufficient to use the non-interacting form 
of the partition function: 



^Z L = ^- 



T ^ 

i 



9i 



Mi 



24vr 2 



Mi / 



7 

+ I5 



d 3 k 
(2^)3 



7TT\ 4 
Mi / 

In (l + e~^~ 



+ In 1 + e 



-/3(e M +/v) 



(6) 



where V is the volume and (3 = T -1 is the inverse temperature. The first term gives the 
contribution of massless particles and antiparticles. Since fj, e 3> m e in all regimes considered 
here, this term applies to both electrons and neutrinos. The degeneracies, are 2 and 1 
for electrons and neutrinos, respectively. The second term gives the muon contributions, 
with degeneracy 



2, = yj k 2 + m 2 and the muon chemical potential designated by 
The pressure, density and energy density of the leptons are obtained from Eq. @ in the 
standard fashion. The total partition function, Z tota \ = ZhZl, where Zh is the partition 
function of the hadrons discussed below. 



4 Stellar Composition and Structure 

4.1 Models of hot and dense hadronic matter 

The properties of neutron stars can be obtained from the well-known hydrostatic equi- 
librium equations of Tolman, Oppenheimer and Volkov [jTj], once the equation of state is 
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specified. At very low densities (n < O.OOl fm~ 3 ), we use the Baym-Pethick-Sutherland flT5| ] 
EOS, while for densities 0.001 < n < 0.08 fm~ 3 , we employ the EOS of Negele and Vau- 
therin |16| . At higher densities the EOS depends on the nature of the strong interactions. 
These are not yet known with certainty, although several intriguing possibilities are cur- 
rently being investigated. In view of this, we will investigate the influence of thermal 
effects and neutrino-trapping on the structure of neutron stars by considering widely dif- 
fering models of dense matter and trying to identify the common features shared by these 
models. 

The models considered include (a) a generalization of a schematic potential model based 
on the work of Prakash, Ainsworth and Lattimer fll7|], which reproduces the results of more 



microscopic calculations [18] of dense matter, (b) a relativistic field theoretical model [19], 
based on the archetypal Walecka model [20], in which baryons interact via the exchange of 
a-, p- and w-mesons, (c) a model based on the chiral Lagrangian of Kaplan and Nelson [^] 
and a model based on a meson exchange picture, in which a kaon condensate occurs at 
about 4no, and finally, (d) a model which allows both quark and hadron phases to be 
present. The composition of the star differs among these models due to differences in the 
nuclear interactions at high density, and, also whether or not strange baryons or strange 
mesons are included in the description of matter. For example, when only nucleons are 
included in models (a) and (b), the proton fraction in matter is determined by the density 
dependence of the symmetry energy; the more rapidly the symmetry energy increases with 
density, the greater is the proton fraction. However, neutrons remain the most abundant 
species in such stars. In model (b), the inclusion of hyperons, which carry strangeness, has 
the effect of substantially softening the EOS at high density. In particular, the presence 
of a substantial number of negatively charged particles, such as the S~ hyperon, raises 
the proton concentration in neutrino free matter and reduces the lepton concentrations. 
In model (c), K~ mesons in the condensate effectively replace electrons to achieve charge 
neutrality, with the result that nearly as many protons as neutrons are found in the dense 
interior regions of the "nucleon" star. In model (d), additional negative charge is provided 
by d and s quarks. 

The composition of matter is significantly altered when neutrinos are trapped. This 
is due to the fact that at the onset of trapping, Yi e ~ 0.4, and the electron fraction Y e is 
significantly larger than that found in a cold catalyzed star. The corresponding changes in 
the structure are quantitatively different among the different models. However, the effects 
of neutrino trapping in matter containing negatively charged, strongly interacting particles, 
either T,~s in hyperonic matter, K~s in Bose condensed matter, or quarks in matter that 
has undergone a phase transition at high density, are qualitatively similar. Since the size 
of thermal effects depends on the relative concentrations, which determine the degree to 
which each constituent is degenerate (T/Tp i <C 1, where Tp i is the Fermi temperature of 
species i), the structural changes are expected to be accordingly different for the different 
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models. 



4.2 Potential models 



Based on a two-body potential fitted to nucleon-nucleon scattering, and a three body term 
whose form is suggested by theory and whose parameters are determined by the binding 
of few body-nuclei and the saturation properties of nuclear matter, Wiringa et al. |jTq| 
have performed microscopic calculations of neutron star matter at zero temperature. At 
high density, there are uncertainties in the three-body interactions, which are reflected 
in the density dependence of the symmetry energy. Calculations at finite temperature 
to encompass an entropy/baryon in the range 5 = 1 — 2 are not yet available. In this 
section, we therefore outline a schematic potential model which is designed to reproduce 

|) of both 



the results of the more microscopic calculations (see, for example, Refs. |18 



nuclear and neutron-rich matter at zero temperature, and which can be extended to finite 



temperature [23]. 

We begin with the energy density 

_ (kin) 



+ 4 Wn > +V(n n ,n p ,T), 



where n n (n p ) is the neutron (proton) density and the total density n 
contributions arising from the kinetic parts are 

d 3 k 



Skin) _|_ (kin) 
~"n. * ^"n 



h 2 k 2 

2m 



(fn + fp) 



(7) 
The 

(8) 



(2tt) 3 

where the factor 2 denotes the spin degeneracy and fa for i = n,p are the usual Fermi- 
Dirac distribution functions and m is the nucleon mass. It is common to employ local 
contact interactions to model the nuclear potential. Such forces lead to power law density- 
dependent terms in V(n). Since repulsive contributions that vary faster than linearly give 
rise to acausal behavior at high densities, care must be taken to screen such repulsive 



interactions [17|. Including the effect of finite-range forces between nucleons, we parametrize 
the potential contribution as 



V(n n ,n p ,T) = \An [§ - (± + s ) (1 - 2x) 



+- 



|-Bn 


3 
2 


- (|+x 3 ) (l-2x) 2 " 


u 


(7+1 


1 + §£'n 


§-(l+x 3 ) (l-2x)2" 





=1,2 



(2Q + 4Zi)2 



d 3 k . 
J2^ 9{ 



1 A i )(/ n + / p 



+{C t - 8Zi) 2 



d 3 k 
(27p 



g{k,Ai)[f n {l-x) + f p x] 



(9) 
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where x = n p /n and u = n/no, with uq denoting equilibrium nuclear matter density. 
The function g(k,Ai) is suitably chosen to simulate finite range effects. The constants 
A, B,a, C\, C2, and B' , which enter in the description of symmetric nuclear matter, and 
the additional constants xq, £3, Z\, and Z2, which determine the properties of asymmetric 
nuclear matter, are treated as parameters that are constrained by empirical knowledge. 

Various limits of the energy density are of interest, and are listed below. Setting x = 1/2 
and f n = f p , the energy density of symmetric nuclear matter is 

„ f d 3 /c Ti 2 k 2 „ , A 2 Bn u a+l 
£nm = 4 J ° U + 1 + B>u°-i 

r cfik 

+ u E C i A 77T^9(k,A t )fn, (10) 

i=l,2 J [ - Z7T) 

and, with x = 0, the corresponding result for pure neutron matter is 

d 3 k Ti 2 k 2 , x . M ,2 §5ra (l - x 3 )u a+1 



n [ d°k h*k z , 1 , M , 2 
7 (2vr) 3 2m J 3 UV u; 1 + §S'(1 - xg)^" 1 

+fu£ " 4 ^) 2 / 7^)3 S(k, K)fn ■ (11) 



i=l,2 



At zero temperature, fi = 9(kF i — k), where kp i is the Fermi momentum of particle i. Thus 
the kinetic energy densities are 

e (^ n ) = |^ 0) n n 5/3 for nuclear matter 

e W = 2 2 / 3 (f^non 5 / 3 ) for neutron matter, (12) 



where = (hkP) 2 /2 m is the Fermi energy of nuclear matter at the equilibrium density. 
To simulate finite range effects, we investigate two commonly used functional forms for 
g(k,Ai), which lead to closed form expressions at zero temperature. 

(i) g(k,Ai) = [1 + (k/Ai) 2 }" 1 : In this case, the finite range termsfl at T = in Eq. @ 
and Eq. (0) may be written as 



u l/3 „l/3' 



V<£> = 3n n C i R i ( V" " arctan ^ ) ( 13 ) 

i=l,2 V ^» 



=1,2 




= S«o« £ (3Q " 4^)^ 3 ( " ^ctan , (14) 



In dynamical situations, such as heavy-ion collisions, it is more appropriate to model matter using 
momentum dependent Yukawa interactions, as pointed out in Ref. j24j. For static matter, both cold and 
hot, the simpler form chosen here is adequate insofar as identical physical properties may be recovered with 
suitable choices of the parameters entering the description of the EOS. 
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where R{ = Ai/(hk F 
in Eq. (£(]) and Eq. 
(u)g(k,Ai) = l- 



Note that the potential energy density for symmetric nuclear matter 



' |jg| ) is the same as that employed in Ref. [17]. 

-(k/Ai) 2 : Here, the finite range interactions are approximated by effec- 
tive local interactions by retaining only the quadratic momentum dependence. The energy 
densities in Eq. (|lG|) and Eq. (11) then take the form of Skyrme's effective interactions [25]. 
The finite range terms, again at T = 0, now read 

3u 2 / 3 



yUr) 



V {fr) 



E °i 



9 2 

^n u 



=1,2 

E 

i=l,2 



1 



(3Q - 4Zi 



1 



3(2«) 2 / 3 
5 



(15) 



Note that, at high density, the quadratic momentum dependence inherent in Skyrme-like 
(SL) interactions will eventually lead to an acausal behavior due to the w 8//3 dependence of 
the energy densities. This situation does not occur in case (i). 

The parameters A, B,a, Ci, C2, and B' , a small parameter introduced to maintain 
causality, are determined from constraints provided by the empirical properties of symmetric 
nuclear matter at the equilibrium density no = 0.16 fm~ 3 . With appropriate choices of 
the parameters, it is possible to parametrically vary the nuclear incompressibility Kq so 
that the dependence on the stiffness of the EOS may be explored. Numerical values of 
the parameters, appropriate for symmetric nuclear matter, are given in Table 1. EOSs 
based on Eqs. (|7[) and (||) are hereafter referred to as the BPAL EOSs if the function 
g(k,Ai) = [1 + {k/Ai) 2 ]- 1 , and the SL EOSs if g(k,A{) = 1 - (k/Ai) 2 . 

In the same vein, by suitably choosing the parameters xo, x$, Z%, and Z2, it is possible 
to obtain different forms for the density dependence of the symmetry energy S(n) defined 
by the relation 



E(n,x) =e(n,x)/n = E(n, 1/2) + 5(n)(l - 2x) 2 + 



(16) 



where E is the energy per particle, and x = n p /n is the proton fraction. Inasmuch as the 
density dependent terms associated with powers higher than (1 — 2x) 2 are generally small, 
even down to x = 0, S(n) adequately describes the properties of asymmetric matter. The 
need to explore different forms of S(n) stems from the uncertain behavior at high density 



and has been amply detailed in earlier publications [17, 2J]. We have chosen to study three 
cases, in which the potential part of the symmetry energy varies approximately as y/u, u, 



and 2u 2 /(l + u), respectively, as was done in Ref. [17|. Numerical values of the parameters 
that generate these functional forms are given in Table 2 for the BPAL and SL EOSs. The 
notation BPALnin2 and SLnin-2 is used to denote different EOSs; n\ refers to different 
values of Kq, and n% = 1, 2 and 3 indicate, respectively, \/u, u and 2ti 2 /(l + u) for the 
dependence of the nuclear symmetry potential energy on the density. 
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The main advantage of casting the schematic parametrization of Ref. [|17]] in the form of 
Eq. (0) through Eq. (||) is that it is now possible to study asymmetric matter at finite tem- 
perature. As an illustration of the calculational procedure at finite temperature, consider 
first the case of pure neutron matter. The evaluation of the baryon density 



n = n r< 



d 3 k 



1 + exp 



i -i 



requires a knowledge of the single particle spectrum 
fi 2 k 2 



2m 



+ U(n,k;T) 



(17) 



(18) 



where the single particle potential U(n,k;T), which is explicitly momentum dependent, is 
obtained by a functional differentiation of the potential energy density in Eq. (|Tl|), with 
respect to the distribution function f n . Explicitly, 



U(n, k; T) = U(n; T) + §u £ (3Q - 4Z<) 



i=l,2 



, k^ 2 

1 + 1 r 



(19) 



for the BPAL EOS, where the explicit momentum dependence is contained in the last term. 
The momentum-independent part is given by 



U(n;T) 



iA(l — Xq)u 



+ 



B(l -x 3 )u a 



1 + §B'(1 - xs)^- 1 



(a + 1) + -B'(l - X3K- 1 



2 , „ \ f d 3 k' fk'\ 

+ ^E(3Q-^)2/^[l + (-) 



U/X2T-2 



fn{k') 



(20) 



For a fixed baryon density n and temperature T, Eq. ( |T7| ) may be solved iteratively for the 
as yet unknown variable 



V 



(21) 



The knowledge of r/ allows the last term in Eq. (20) to be evaluated, yielding U, which may 
then be used to infer the chemical potential from 



Hn=Tr]-U . 



(22) 
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which is required in the calculation of the single particle spectrum in Eq. (jl8|). With 
this ek, the energy density in Eq. ([ll]) is readily evaluated. The entropy density has the 
same functional form as that of a non-interacting system: 

s = ~2 J ^3 [/„ln / n +(l-/n)ln (1 - /„)] , (23) 

from which the pressure is obtained using 

p = s T + nfi n -e. (24) 

The above procedure is also applicable, with obvious modifications, to a system contain- 
ing unequal numbers of neutrons and protons, which is generally the case for charge-neutral 
matter in beta equilibrium. 



4.2.1 Pure neutron matter 

The influence of finite entropy on the structure of a star is most easily studied when the 
star is idealized to be composed of neutrons only. The top panels in Fig. 2 show the density 
dependence of the Landau effective mass, m* = kp/{de^Jdk)\h F , for the BPAL22 and SL22 
EOSs. For the BPAL22 model considered, the effective mass is weakly dependent on density 
and has a value ~ 0.6m, where m is the bare mass. On the other hand, for the SL22 EOS 
the variation of the effective mass with density (top panels) is significantly more rapid. 
This qualitative difference follows from the explicit expressions for the effective mass: 
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for the BPAL EOS, and 



(25) 
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(26) 



for the SL EOS. 

Turning to the pressure (center panels of Fig. 2), it is possible to analytically estab- 
lish the quadratic dependence on the entropy per baryon, and, perhaps more importantly, 
understand the order of magnitude of the increase in the maximum mass if we analyze the 
thermal contributions using the methods of Fermi liquid theory. Following the analysis in 
Ref. |^]], the thermal pressure of interacting neutrons may be cast in the form 
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where the Fermi temperature, Tp = h 2 k F /2m*, sets the scale of the temperature depen- 
dence of the thermodynamical functions. In the dense central regions of the star, T/Tp -C 1, 
so the neutrons, which are the only constituents considered now, are in a highly degenerate 
configuration. Since, in the degenerate limit, the entropy per particle S = (tt 2 /2)(T/Tp), 
the ratio of the thermal pressure to that of the zero temperature pressure Po niay be ex- 
pressed as 



Pth 
P 



o 



— S" 2 + 



1 



3 d In m* 



2 dlnn 



j _|_ Ppot 



"I -1 



Phi 



kin 



(28) 



where Pkin and P po t, the kinetic and potential pressures, as well as the Landau effective 
mass, m* , in this relation refer to zero temperature matter. This equation adequately 
reproduces the exact thermal pressure at the entropies likely to be relevant in the evolution 
of a neutron star. The differences of Eq. (^) from the ideal gas result are obvious. First, 
the density dependence of the effective mass introduces a correction; this is clearly the 
origin of the larger thermal pressure for the SL22 case compared to the BPAL22 case for a 
given entropy (see Fig. 2). Second, the pressure arising from potential interactions, which is 
generally larger than the kinetic pressure due to the predominance of repulsive interactions 
at high density, produces a significant reduction in the ratio P t h/Po- 

Since the thermal pressure increases quadratically with the entropy/baryon, S = s/n, 
the neutron star masses should show a similar behavior, as we see from the bottom panels 
of Fig. 2. Qualitatively similar results are obtained for our other EOSs with different 
behavior for the symmetry energy. Quantitative results for the basic physical attributes 
of a maximum-mass star are shown in Table 3 for the EOSs termed BPAL22 and SL22, 
respectively. It is clear that the maximum mass at finite entropy is well approximated by 

M max (S) = M max (0) [l + \S 2 H J , (29) 

where the coefficient A is EOS dependent. The values of A given in Table 3 are quite small, 

~ icr 2 . 

These results highlight the point that the magnitude of the increase in the maximum 
mass is chiefly governed by the magnitude of the thermal pressures at fixed entropy. In the 
simple models considered above, the thermal pressures depend sensitively on the momen- 
tum dependence of the nuclear interactions. Insofar as one can establish this momentum 
dependence, either from experiment or theory, the thermal pressures may be constrained. 
For example, the energy dependence of the real part of the optical model potential required 
to explain proton-nucleus scattering experiments provides a stringent constraint on the 
momentum dependence of the single particle potential at nuclear density, but for values of 
x ~ 0.5. For higher densities in the range (2 — 3)no attained in GeV/particle heavy-ion 
collisions, the momentum dependence largely governs the flow of matter, momentum and 
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energy p4p . While the BPAL single particle potentials are consistent with the observed 
behavior, the SL potentials, with their quadratic momentum dependence, are known to 
be inconsistent with data at high momentum |[28|| . Also, at high density, the quadratic 
momentum dependence has the disadvantage that it leads to an acausal behavior, due to 
the u 8 / 3 dependence of the energy densities in Eq. (|T5|), 



4.2.2 Neutrino-free matter in beta equilibrium 



We turn now to the more realistic case in which matter consists of neutrons, protons, 
electrons, and muons, with their relative concentrations determined from the conditions of 
charge neutrality, Eq. (pi), and equilibrium under beta decay processes in the absence of 
neutrino trapping, Eq. (|3|). The EOS of strongly interacting matter above nuclear density 
is given by Eq. (Q) and Eq. (||). The EOS of leptons is obtained from Eq. 

The relative concentrations, the electron chemical potential, and the individual en- 
tropies per baryon (for S = 1) are shown in Fig. 3 for the BPAL22 and SL22 EOSs, both 
of which are characterized by a symmetry energy whose potential part varies linearly with 
density. Note that at high density, the proton concentration lies in the range (20 — 30)% 
(for these EOSs), which is balanced by an equal amount of negatively charged leptons to 
maintain charge neutrality. Also, the lepton contribution to the total entropy per baryon 
is comparable to that of the degenerate nucleons at high density, thus lowering the nucleon 
contribution at fixed S. This results in a smaller increase in the maximum mass than when 
leptons are not present, as is evident from Fig. 4, where the nucleon effective masses, isen- 
tropic pressures, and the mass curves for the corresponding EOSs are shown. The increase 
in the maximum mass is again quadratic with entropy. This is due to the fact that both 
nucleonic and leptonic pressures stem mostly from quadratic terms in the entropy, terms 
involving higher powers of entropy giving rather small contributions. When the nucleons 



are degenerate, the generalization of Eq. (28) is 
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(30) 



where Y{ = ni/n and Tp t = (ft/cpj 2 / '(2m*) denote the number concentration and the Fermi 
temperature of species i, respectively. Further, Pkin, Ppot, and Po = Pkin + Ppot are the 
kinetic, potential, and the total pressures at zero temperature and include the contributions 
from both neutrons and protons. 

Quantitative results for the physical attributes of the maximum mass star are summa- 
rized in Table 4 for the different BPAL EOSs which have varying stiffnesses and different 
density dependence of the symmetry energy, the latter determining the proton and lepton 
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concentrations in the star. Results for the SL interactions with a linear potential contri- 
bution to the symmetry energy, but for EOSs with different stiffness are shown in Table 5. 
The results in these tables reflect the influence of entropy on the gross properties of stars. 
With few exceptions, the temperature at the core remains below 100 MeV. Although the 
central density and pressure of the maximum mass star are significantly reduced, A is uni- 
formly ~ 10~ 2 , so the increase in the maximum mass amounts to only a few percent of the 
cold catalyzed star, even up to S = 2. In all cases studied, the increase in mass is quadratic 
with entropy. 

The moments of inertia, /, displayed in Tables 4 and 5 refer to maximum mass stars, 
so that the effect of cooling a given star from entropy per baryon, S, of 2 to cannot 
be directly assessed. We therefore show in Fig. 5 the moment of inertia as a function of 
the central density and of the baryonic mass, Mb- (The full dots on the various curves 
represent the maximum mass configurations.) Since the baryonic mass is proportional to 
the number of nucleons, it will remain fixed as the star cools, in the absence of accretion. 
Fig. 5 shows that the moment of inertia decreases as the star cools from S = 2 to S = 0. 
The magnitude is dependent on the mass in question, but a typical value is ~ 20%. Since 
angular momentum is conserved, again in the absence of accretion, this implies an increase 
of the angular velocity of similar order. Thus a significant spin-up of the neutron star is 
expected to occur during the cooling process. 

4.2.3 Neutrino-trapped matter in beta equilibrium 

We turn now to the case in which neutrinos are trapped in matter. As mentioned earlier, we 
fix the electron lepton number at Yj, e = 0.4 and the muon lepton number at Yl^ = 0. Fig. 6 
shows the effects of neutrino trapping on the relative abundance, the chemical potentials, 
and the partial entropies for an entropy per baryon S = 1. The major effect of trapping is 
to keep the electron concentration high so that matter is more proton rich in comparison to 
the case in which neutrinos are not trapped (cf. Fig. 3). The EOS with trapped neutrinos 
is less stiff than that without neutrinos, since the decrease in pressure due to the nuclear 
symmetry energy is greater than the extra leptonic pressure. This is reflected in the limiting 
masses obtained with neutrino trapped matter (see Table 6), which are lower than those 
without neutrinos (cf. Table 4). Results for the maximum mass stars are summarized in 
Table 6 for the BPAL models with different stiffnesses. The decrease in the maximum mass 
due to the effect of trapping is generally larger than the increase due to thermal effects for 
S = 1. For S = 2 the thermal contributions are larger, which results in a delicate balance 
between the two competing effects. 

It is instructive to contrast the effects of thermal contributions on white dwarf and 
neutron star structures. Unlike the case of a white dwarf, it is not possible to analytically 
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predict the effects of finite entropy on the structure of a neutron star. While both configu- 
rations are degenerate, and thus one expects the dominant effects of finite entropy to enter 
quadratically, the role of interactions in the two cases are quite different. In contrast to the 
case of a white dwarf, in which the equation of state is highly ideal [jDJ, the equation of 
state in a neutron star is strongly non-ideal. This results in a white dwarf having a much 
greater sensitivity to entropy than a neutron star. For a white dwarf, in which the electron 
pressure dominates, the thermal correction to the Chandrasekhar mass is about 10% at an 
entropy per baryon of 1 |29|. In neutron stars, the pressure support is largely provided by 
the strongly interacting baryons, which have relatively smaller thermal contributions to the 
pressure and, therefore, a smaller increase in the maximum neutron star mass. As a result, 
the compositional variables of the EOS play a more important role than the temperature 
for the structure of neutron stars. 



4.3 Field theoretical models 

From the results of the previous section, we notice that the central density of the maximum 
mass star typically exceeds (4 — 5) hq. At such densities, the Fermi momentum and effective 
nucleon mass are both expected to be on the order of 500 MeV. Thus a relativistic descrip- 
tion is preferred. Relativistic local quantum field theory models (see, for example, Ref. J2(J) 
of finite nuclei and infinite nuclear, and neutron star, matter have had some success, albeit 
with rather more schematic interactions and with less sophisticated approximations than 



their nonrelativistic counterparts (see, for example, Refs. [18, ^]). It is our purpose here 
to examine the effects of finite entropy and neutrino-trapping in a relativistic description 
and to contrast the results with those of potential models. 

Specifically, we employ a relativistic field theory model in which baryons, B, interact 
via the exchange of a-, p-, and w-mesons. In the case that only nucleons are considered, 
B = re, p; this is the well-known Walecka model ^0|, which we evaluate in the Hartree 
approximation (or, equivalently, at the one- loop level). It has been shown, however, that 
hyperons significantly soften the zero-temperature equation of state p0| , |3l| . Therefore, we 
shall also consider the case where the hyperons, A, E, and S, are included in the set of 
baryons B. (The inclusion of the spin-| A quartet and the Q~ hyperon is straightforward, 
but does not quantitatively alter the results, since they appear at densities much higher 
than found in the cores of stars.) Specifically, our Lagrangian is 

C H = Y,B (i^d^ - g^Bl^n - 9 P B^h^ ■ t - M B + g aB a) B 

B 

+\d lx ad> l cr - \mla 2 - U{a) 
-\F^ + \ml^ 

-iB^-B^ + im^-b^. (31) 
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Here the p- meson field is denoted by b^, the quantity t denotes the isospin operator which 
acts on the baryons, and the field strength tensors for the vector mesons are given by the 
usual expressions:- = d^uj u - d u uj^, = d fJi b u - d u h^. 

It is straightforward to obtain the partition function for the hadronic degrees of freedom, 
denoted by Zjj, 



hxZ H = pV 



\ml<4 + \m 2 p bl - \rn\a 2 - U{a) - £ AE(M* B ) 



B 



+^E/(03^(l + e-W--)) . (32) 

We shall consider the relativistic Hartree approximation, for which we take U(a) = 0. 
We shall also consider the mean field approximation, in which the negative energy sea 
contributions AE(M B ) is neglected and additional scalar self-couplings are included with 
U(a) = (6M/3)(5 -at(t) 3 + (c/4)(<7 <J ivO') 4 - The contribution of antibaryons are not significant 
for the thermodynamics of interest here, and is therefore not included in Eq. (|32j). Here, the 



effective baryon masses M B = M B — g a B& and E* B = yk 2 + M* B . The chemical potentials 
are given by 

PB = »B+ QuiB^O + QpBhBbo , (33) 

where t%B is the third component of isospin for the baryon. Note that particles with 
t-3B = 0, such as the A and E° do not couple to the p. When hyperons are included, the 
negative energy sea contribution from all baryons, inclusive of the hyperons, is considered 
as indicated by the notation AE(M B ). 

The shift in the energy density of the negative energy baryon states is evaluated in the 
one loop Hartree approximation. After removing divergences, AE(M B ) can be written |Q| 
in the form 

AE(M* B ) = - JL Ul - £ + In £) M B (M B - M B f - In £ (M B - M B ) A 



M M J °' M 

M* 

+M* B 4 In -/ + M 3 B (M B - M* B ) - ImUMb - M* B f 

M B 

+fM B (M B - M* B f - §(M B - M B f 



(34) 



where M = 939 MeV is the nucleon mass. Here, the necessary renormalization introduces 
a scale parameter, p r . For the standard choice [20|, 33] of p r /M=l (termed RHA), the first 



two terms in Eq. (|34| ) vanish. This will not be the case for other choices of p r /M, which 
introduce explicit cr 3 and a 4 contributions. At the phenomenological level, the a 3 and 
cr 4 couplings, generated from the baryon-loop diagrams, modify the density dependence of 
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the energy, which makes it possible [32] to obtain nuclear matter compression moduli that 
are significantly lower than in the standard RHA. We shall exploit this freedom to vary 
H r /M and we call this approach the modified relativistic Hartree approximation (MRHA). 
In previous work ||l9|| without hyperons, we found that while neutron star masses do not 
significantly constrain fj, r /M, finite nuclei favor a value of 0.79. However, our interest here 
is to explore the dependence on the varying stiffness of the equation of state, which is 
permitted by modest variations in fj, r /M, for the purpose of studying the impact on the 
structure of the star at finite entropy. 

Using Zh, the thermodynamic quantities can be obtained in the standard way. The 
pressure Ph = TV^ 1 In Zh, the number density for species B, and the energy density ejj 
are given by 

d 3 k / .A- 1 



e H = \rnla 1 + U(a) + i m £wg + 1^2 + &E{M%) 
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^[^^(e^i-^ + l)- 1 . (35) 
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(27T) 

The entropy density is then given by sh = P(£h + Ph ~ J2b P^b^b)- 

The meson fields are obtained by extremization of the partition function, which yields 
the equations 

m S w o = ^2 9wB n B ; m 2 p b = ^2 gpBhsnB , 

B B 

m> .jm + ^j^ ^ + r+ _jl [A£{MM } . m 

The additional conditions needed to obtain a solution are provided by the charge neu- 
trality requirement, Eq. (^), and, when neutrinos are not trapped, the set of chemical 
potential relations provided by Eq. (||). For example, when i = e~ , this implies the 
equalities 

MA = Ms = Ms = Mn ! 
Me- = Ms- = Mn + Me , 

Mp = ME+ = Mn - Me • (37) 
In the case that the neutrinos are trapped, Eq. @ is replaced by Eq. (||). The new 
equalities are then obtained by the replacement fi e — ► \x e — fi Ue in the above equations. The 
introduction of additional variables, the neutrino chemical potentials, requires additional 
constraints, which we supply by fixing the lepton fractions, Yu, as noted above. 
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4.3.1 Neutrino- free, non-strange baryonic matter 

For the MRHA calculations, the nucleon coupling constants were fitted to the equilibrium 
nuclear matter properties: a binding energy per particle of 16 MeV, an equilibrium density 
of 0.16 fm~ 3 , and a symmetry energy of 30 MeV (at the equilibrium density), as in Ref. |J. 
Defining Cf = (guv M / 'm;) 2 , where N represents n or p, the parameters C 2 , C 2 , and C 2 p are 
reproduced in Table 7 for a range of values of \i r jM. We see that they encompass a fairly 
wide range of values for the nucleon effective mass, Mjy-, and the compression modulus, 
Kq, at saturation. Correspondingly, the stiffness of the equation of state shows significant 
variations. Note that Kq alone is not always a reliable indicator of the stiffness at high 
density; for example, the equation of state for /j, r /M=0.73 is noticeably stiffer than for 
fi r /M=1.25, yet the compression modulus is smaller, see Ref. [jOJ. 

In the mean field approximation, we use one of the parameter sets of Glendenning and 
Moszkowski |}4|], so that comparisons may be made. Specifically, we use the couplings 

b = 0.008659 c = -0.002421 
Cl = 109.14, C 2 a = 224.78 and C 2 p = 108.49 . (38) 

With these constants, the equilibrium density of nuclear matter is no = 0.153 fm~ 3 with 
the Landau effective mass M^ sat /M = 0.827 and compression modulus Kq = 240 MeV. 
This model is referred to as the GM model. 

Table 8 shows the basic properties of stars in beta equilibrium for the MRHA model 
with various values of fx r /M and for the GM model. The structural changes at finite entropy 
compared to the zero entropy case are qualitatively similar to those given by the BPAL 
potential model in the previous section. In most cases, A is a little smaller in magnitude, 
and this is probably due to the strong repulsion of the w-meson at high density, which 
increases the P po t term in Eq. ([28|). Again, a nearly quadratic increase of the maximum 
mass with entropy is obtained. The central temperatures, T c , are lower here than in the 
non-relativistic case. Since it is the entropy per baryon that is constant, the temperature 
will vary with density. The temperature profile as a function of density is shown by the full 
curves in the upper panel of Fig. 7 for /i r /M=1.25. The temperature is a maximum at the 
center of the star (here the density ratio is ~ 7 for a maximum mass star, see Table 8) and 
decreases with decreasing density, the fall off being particularly marked at low density. The 
density of a neutron star is approximately constant in the interior and drops to zero over a 
radial distance AR/R ~ 0.1. Thus the interior of the star will have a constant temperature, 
and this will fall off rapidly in the surface region. 

The top panel in Fig. 8 shows the composition of the star for the model with ^ r /M=1.25. 
The middle panel shows the electron (or muon) chemical potential, while the bottom panel 
shows the individual contributions to the total entropy. Note that the lepton contribution 
to the entropy is smaller than in the non-relativistic case. In Fig. 9, we show in the top panel 
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the behavior of the Landau effective masses for neutrons and protons; these are defined by 

37r 2 n r 



m 



E 



Fn 



+ M* 2 , where k Fn 



for neutrons, and similarly for protons. 



The isentropic pressures are shown in the middle panel of Fig. 9, and the bottom panel 
indicates the star mass versus central density ratio for fixed entropies. These results are 
quite similar to those shown in Figs. 3 and 4 for BPAL22, although the maximum masses 
obtained here are somewhat larger. 

As with the potential models, the nucleonic contributions to the thermal pressure may 
again be simply estimated in the degenerate limit, T/Tp <C 1. For one nucleon species, the 
Fermi temperature Tp = k F /(2E F ), and the thermal pressure for this relativistic model is 
given by an expression of the form (see Ref. \27\ for more details) 
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(39) 



where the term containing the logarithmic derivative arises from the density dependence 
of the effective mass. The entropy per particle in the degenerate limit is again given by 
5 = (it 2 /2)(T/Tp), but now with the Fermi temperature appropriate for a relativistic 
spectrum. Generalizing to two nucleon species, the total thermal pressure of the nucleons, 
up to quadratic terms in 5, may be written as 
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(40) 



where = n^/n with n the total nucleon density. This provides an accurate approximation 
to the exact results. 



4.3.2 Neutrino-trapped, non-strange matter 

With Yp e = 0.4 and Ylu = 0, Fig. 10 shows the effects of neutrino trapping on the relative 
abundance, the chemical potentials, and the partial entropies for an entropy per baryon 
5 = 1. The temperature profile of the star, indicated by the dotted curves in Fig. 7, differs 
rather little from the untrapped case. Results for the maximum mass stars are summarized 
in Table 9 for the MRHA and GM models. Note that the effect of trapping is similar to that 
found earlier with the potential model, especially for the abundance of protons due to the 
large concentration of electrons. The amount by which the maximum decreases, ~ 0.07Mq 
when neutrinos are trapped, is also similar to the results of the potential models. Since 
thermal effects are smaller in the relativistic models, the maximum mass of the 5 = 2, 
neutrino-trapped star is always less than that of the cold 5 = 0, neutrino-free star. 
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Together with the findings of the earlier section, we conclude that in matter in which 
the only baryons are neutrons and protons, neutrino trapping usually decreases the maxi- 
mum mass by a larger amount than thermal effects increase it. 



4.3.3 Neutrino-free, strangeness-rich baryonic matter 

When hyperons are included, their coupling constants are needed; however, these are largely 
unknown. We shall assume that all the hyperon coupling constants are the same as those of 
the A, for which we can take some guidance from hypernuclei. Following Glendenning and 



Moszkwoski [34], the binding energy of the lowest A level in nuclear matter at saturation 
yields 

B\ = fiA- M A = x w &„w + Ml - M A , or 
-28 MeV = x^g^nUo - x a g an a , (41) 

where x a = g a Aldan arid x w = g^A/gum- Further, as suggested in Ref. [34] on the basis of 



fits to hypernuclear levels and neutron star properties, we take x a = 0.6. The value of x^ 
may then be determined from Eq. (|4l|). For the p meson, we take x p = g p A/9pn = x a . The 
alternative choice, x p =i w , is found to produce essentially similar results. 

In the MRHA model, the negative energy sea contributions AE(M^) from all baryons, 
inclusive of the hyperons, contribute even when the positive energy states of the hyperons 
are empty. This entails a redetermination of the constants for the case in which hyperons 
are included. In Table 10, the coupling constants that reproduce nuclear matter saturation 
properties and the binding energy of the lowest A level in nuclear matter at saturation 
(with x a = 0.6) are given for different choices of [i r jM. Also shown are the nucleon 
effective masses and the compression moduli. Since the negative energy sea contributions 
of all the baryons are positive, the values of are somewhat smaller than those in Table 
7 for the case in which only nucleons are considered. Consequently, the values of Kq are 
also smaller than those in Table 7. 

For the MRHA model, with // r /M=1.25, we show in Fig. 11 the relative fractions (top 
panel) of the baryons and leptons in beta equilibrium, and the electron chemical potential 
(bottom panel) as a function of baryon density at zero temperature. One expects that A, 
with a mass of 1116 MeV, and the £~, with a mass of 1197 MeV, first appear at roughly 
the same density, because the somewhat higher mass of the S~ is compensated by the 
the presence of the electron chemical potential in the equilibrium condition (see Eq. (|37|)) 
of the £~. More massive, and more positively charged, particles than these appear at 
higher densities. Notice that with the appearance of the negatively charged X~ hyperon, 
which competes with the leptons in maintaining charge neutrality, the lepton concentrations 
begin to fall. This is also reflected, for example, in the magnitude of the electron chemical 
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potential which saturates at around 200 MeV and begins to fall once the S~ population 
begins to rise rapidly. (That the negatively charged S~ is the cause for fi e to fall with 
density may be verified by allowing only the neutral particles to appear. In this case, \i e 
continues to rise with density, albeit with a reduced slope compared to the case in which 
no neutral hyperons are present.) The rapid build up of the other hyperons with increasing 
density has two major consequences. First, the system is strangeness-rich at high density, 
with nearly as many protons as neutrons. Second, since at a given total baryon density 
the system contains many more baryon species with sizeable concentrations, the EOS is 
considerably softer than when no hyperons are present. This causes the maximum mass to 
be reduced |3(], 31 1. 

Fig. 12 contains the corresponding results at an entropy per baryon S = 1. As expected, 
at finite temperature the hyperons attain significant fractions at lower baryon densities 
than at zero temperature. The order of appearance of the various hyperons follows from 
the chemical potential relations, Eq. (|37|), and their differing masses. The bottom panel 
of Fig. 12 shows that the baryons carry most of the entropy, since the lepton populations 
remain low at high density, due to the magnitude of the electron chemical potential. 

Table 11 summarizes the gross features of the maximum mass star populated with hy- 
perons. Compared to the case in which only nucleons are present, the addition of hyperons 
causes the central temperatures to be reduced. This is evident by comparison of the bot- 
tom and the top panels in Fig. 7. This figure also shows that with hyperons present the 
temperature changes rather little with density until u < 2, so that a constant temperature 
would be achieved over much of the star. The softening introduced in the MRHA EOS by 
hyperons is evident in the maximum masses in Table 11. These are about 0.4 — 0.9Mq, 
smaller than the results of Table 8, for which only nucleons are allowed in matter. Notice 
that in some cases, the maximum mass falls below 1.44M , and the pressure support of 
finite entropy is not adequate to raise the maximum mass above 1.44M . We also give 
here results for the mean field GM model |34j (with x a = x p = 0.6 and = 0.659). This 
model shows a similar reduction in the maximum mass of ~ 0.5 M@, due to the softening 
induced by hyperons. Qualitatively similar results, albeit with somewhat larger limiting 
masses, are obtained for the other choices of couplings in Ref. p4|]. 

In contrast to the maximum masses of stars containing nucleons and leptons only, the 
maximum masses of the hyperon populated stars do not show a regular behavior with in- 
creasing entropy. In those cases for which the maximum mass increases with entropy, a 
quadratic increase is observed. However, with fi r /M=l, a decreasing trend with entropy is 
found. We have verified that this surprising behavior does not violate any laws of thermo- 
dynamics. It can be traced back to the softening of the EOS caused by the appearance of 
hyperons at relatively smaller densities than are found at zero temperature. 
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4.3.4 Neutrino-trapped, strangeness-rich matter 

Fig. 13 is the counterpart of Fig. 11 for the case in which neutrinos are trapped at zero 
temperature. Trapped neutrinos have a large influence on the charged hyperon thresholds. 
For example, the appearance of the S~ hyperon is governed by /i s - = fi n + fx. Here, 
/i = fi e —[i Ue is much smaller than in the untrapped case for which /x = /i e , so the appearance 
of the S~ occurs at a higher density. With the appearance of hyperons, the neutrino 
population begins to increase with density, in contrast to the monotonic decrease exhibited 
in the hyperon free case. The finite chemical potential of the neutrinos requires the electron 
chemical potential to be at a higher value than in the neutrino free case to maintain chemical 
equilibrium. Thus, unlike the neutrino free case, the electron chemical potential increases 
with density. Notice that muons play little role here. The preponderance of negatively 
charged particles, both leptons and baryons, now has the consequence that the system is 
proton rich over an extended region of density. These qualitative features are retained also 
at finite entropy with baryons carrying most of the entropy, as shown in Fig. 14. 

The physical properties of the maximum mass stars with trapped neutrinos are listed 
in Table 12. The changes due to entropy alone are small and not always in the direction of 
increasing the maximum mass. Notice, however, that for each entropy shown, the maximum 
masses are all about 0.2M@ larger than those in Table 11 for neutrino free stars. Since 
the star has to cool down from an S ~ 1 configuration, with neutrinos trapped, to a 
configuration of S ~ without neutrinos, the maximum stable mass decreases. This may 
be contrasted with the case in which no hyperons are present, in which case neutrino 
trapping and finite entropy effects are opposed to each other and effectively cancel. This 
is a general result that stems from the softening induced by the presence of negatively 
charged hadrons in the star, as demonstrated in the next section, where kaon condensation 
is allowed to occur in dense matter. 

4.3.5 Metastability of neutron stars with strangeness-rich baryons 

The most striking conclusion of the above discussion is the possibility of metastable neutron 
stars if matter contains hyperons. Metastable stars occur within a range of masses near the 
maximum mass of the initial configuration and remain stable only for several seconds after 
formation. In contrast, in matter with only nucleons, any star that is below the mass limit 
of the initial configuration will be stable during the subsequent evolution (in the absence 
of mass accretion). 

In an idealized picture, the two features that govern the evolution of the maximum 
mass are the lepton content and finite entropy. The latter plays a minor role, but effects 
of the high lepton content are very important. The combined effects of lepton fraction 
and finite entropy are shown in Fig. 15. Here, the abscissa is the baryonic mass which is 
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proportional to the number of baryons in the system and is constant during the evolution 
of the star (in the absence of accretion of matter, again, an idealization). The ordinate is 
the gravitational mass, hitherto referred to simply as the mass, which includes interactions 
and, therefore, changes as the star evolves. If hyperons are present (lines ending with a 
dot), then deleptonization, that is the transition from Y^ e = 0.4 to Y v = 0, accompanied by 
heating from S = 1 to S = 2, lowers the range of baryonic masses that can be supported 
by the equation of state (of model GM here) from about 1.95M Q to about 1.73M0. The 
window in the baryonic mass in which neutron stars are metastable is thus about O.22M 
wide. On the other hand, if hyperons are absent (lines ending with a star), the baryonic 
mass increases during deleptonization, and no metastability occurs. Similar results apply 
to the gravitational masses, which can be obtained from the baryonic masses with the help 
of the right panel. 

Fig. 16 illustrates the evolution of the maximum mass during deleptonization. For 
clarity, the small effects of finite temperature are not considered here. When neutrinos 
diffuse out, the neutrino fraction decreases from an initial maximum value of about Y Ue = 
0.08, corresponding to Y^ e = 0.4, to zero. The evolution of the neutrino fraction with time 
can only be obtained by solving the neutrino and energy transport equations, as was done, 
for example, by Burrows and Lattimer |fflj. This is beyond the scope of the present work. 
However, if one assumes that the neutrino fraction decreases approximately uniformly with 
time, then one obtains a good picture of how the maximum mass will evolve. In Fig. 16, 
the maximum mass is shown as a function of the electron neutrino number per baryon, Y Ue , 
for the GM model. It is evident that different compositions lead to different trends. In 
nucleonic matter, the maximum mass slowly increases as neutrinos diffuse out; in hyperonic 
matter, maximum mass decreases. It is also clear that the rate of change is largest shortly 
before the neutrino fraction drops to zero. Similar qualitative behavior is obtained for other 



choices of mean field parameters from Ref. [34| 



4.4 Matter with kaon condensation 

The idea that, above some critical density, the ground state of baryonic matter might 
contain a Bose-Einstein condensate of kaons is due to Kaplan and Nelson [21]. The formu- 



lation, in terms of chiral perturbation theory, was subsequently discussed by Politzer and 
Wise |35|[ and Brown et al. J36|. The composition and structure of kaon condensed stars 
and also some evolutionary aspects were considered by Thorsson, Prakash and Lattimer ||. 
Further related calculations at the mean-field level may be found in Refs. ||, ||. Most 



recently, loop contributions have also been investigated [^9[. Depending on the parameters 
employed, particularly the strangeness content of the proton, kaon condensation is typically 
found to occur at about 4 times the equilibrium nuclear matter density. This has the effect 
of softening the equation of state and lowering the maximum mass of the neutron star. 
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In chiral models the kaon-nucleon interaction occurs directly via four point vertices; 
however, it can also be modelled as an indirect interaction that arises from the exchange of 
mesons [[EJ. The latter approach has the virtue that it is more consistent with the meson 
exchange picture that is usually employed for the baryon interactions, and it is of interest 
to compare the predictions with those of the chiral model. A further question is the role 
that hyperons might play in addition to kaons, since we have seen that they appear at a 
density similar to, or somewhat lower than, the condensate threshold. This raises the issue 
of the interplay between the strange baryons and mesons and the net strangeness content 
of neutron stars. These questions have been the subject of recent work [41, 42, fSJ, and we 
shall address them in the present context. 

Our purpose here is to determine the impact of a kaon condensate on protoneutron 
stars. In order to keep the number of calculations within bounds, we will focus on the 
mean field case with the GM parameters. We will study the effect of neutrino trapping 
in chiral and meson exchange models, both with and without hyperons. As regards fi- 
nite entropy effects, we have seen that they are less significant than neutrino trapping. It 
is not clear how to develop a consistent finite temperature formalism for the chiral case, 
while it is reasonably straightforward for the meson exchange model. We therefore study 
the effects of finite entropy in the latter approach and restrict our chiral formalism to T = 0. 



4.4.1 Chiral formalism 



The Kaplan-Nelson SU(3) x SU(3) chiral Lagrangian for the kaons and the s-wave kaon- 
baryon interactions takes the form 

C K = IfTvd^Ud^U+CTimgiU + U j -2) + iTvB^[V^,B] + ai TrB(£m q £ +h.c.)B 



+a 2 TrBB(£m q £ + /i.e.) + a 3 {TrBB} {Ti(m q U + h.c.)} 



(42) 



Here, U = £ 2 is the non-linear field involving the pseudoscalar meson octet, from which we 
retain only the contributions- 



U 



exp 



'V2i 



M 



M 




(43) 



The baryon octet - nucleons plus hyperons - is given by 



/JIS° + V /±A 



B 



A 



V 
n 



(44) 
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In Eq. (p2y, the quark mass matrix m q = diag(0, 0, m s ); i.e., only the mass of the strange 
quark is taken to be non-zero. For the mesonic vector current, Va, only the time component 
survives in an infinite system with Vq = ^(tfdo^ + £(?o£ )■ Also, the pion decay constant 
/ = 93 MeV, and C,a±,a2 and 03 are constants. After some algebra, the relevant part of 
Ck takes the form 

s ' m X\ 2 L rs+nurs- , * (K+d K~ -K-d K^ 



X J \ 4/ 2 cos 2 ±x 



B 

1 BB 



~[ m K + TH^2\( a ^ + fl 2)( 1 + Vb?b) + (01 - a 2 )(<?B - Yb) + 4a 3 
V Z J b 

+ a 2)(2AA + V3[E°A + AS ])) , (45) 

0/ / COS 2 X J 

where gs and Yb are the baryon charge and hypercharge, respectively, and q B = ^Ys + tss- 
(The hypercharge is given in terms of the baryon number and strangeness by Y B = bs+Ss-) 
In Eq. (^), we have defined x 2 = 2K + K~ / f 2 and taken the kaon mass to be given by 
m 2 K = 2Cm s /f 2 . We have not included in Eq. (]45| ) terms which simply give a constant 
shift to the baryon masses; they indicate that aim s = —67 MeV and a^vtis = 134 MeV, 
using the hyperon-nucleon mass differences. The remaining constant a^m s is not accurately 
known, and we shall use values in the range —134 to —310 MeV corresponding to to 20% 
strangeness content for the proton. The corresponding range for the kaon-nucleon sigma 
term, 

\{a x + 2a 2 +4a 3 )m s , (46) 



2 



is 167-520 MeV. Some guidance is provided by recent lattice gauge simulations [44], which 
find that the strange quark condensate in the nucleon is large, i.e., (N\ss\N) = 1.16 ±0.54. 
From the relation m s (ss) p = — 2(a 2 + a^)m s and using m s = 150 MeV, we obtain a^m s = 
— (220 ± 40) MeV, which is in the middle of our range of values. 

The pure kaon part of Eq. ( f45|) gives a contribution to the grand potential, fix = 
—TIiiZk, which can be evaluated in the mean field approximation by writing the time 
dependence of the fields K± = -^fOe^ §. Here, 6 gives the condensate amplitude, and 
H is the kaon chemical potential. One easily finds 

n K = Vf 2 (2m 2 K sin 2 ±# -\^ 2 sin 2 6<) . (47) 

The kaon-baryon interactions in Eq. (fiq) can be incorporated in the baryon Lagrangian 
of Eq. ( |3l"D by suitable definitions of the effective masses and chemical potentials. Notice 
first that the A — E° mass matrix needs to be diagonalized. We write it in the form 



2a 2(3 

W 27 



where 
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2a = Ma - g a \a + (|ai + |a 2 + 4a 3 )m s sin 2 i( 



2/3 = 3 a(ai + a 2 )m s sin 2 i<9 
Then it is straightforward to obtain the eigenstates and the masses 



2 7 = M s - g alZ a + ( ai + a 2 + 4a 3 )m s sin 2 ±0 . (48) 



H x = — ; H 2 = iy ^ UZj W ith 2(35 = 7 - a - J( 7 - a) 2 + 4/? 2 

(l + <5 2 )5 (l + <5 2 ): V 



)2 



= 7 + a + v /( 7 -a) 2 + 4/? 2 ; M^ 2 = 7 + a - ^( 7 -a) 2 + 4/3 2 . (49) 

Since 7 > a here, in the limit of no mixing (/3 = 0) 6 is zero. Henceforth, the sum over 
baryon states, B, includes H\ and H 2 along with the n, p, H 0, ~. For the latter cases, 

the masses are given by 

M* B = M B - g aB a + [{a x + a 2 )(l + Yb9b) + (ai - a 2 )(g B - y B ) + 4a 3 ]m s sin 2 ±6> .(50) 

The chemical potentials //b are given in terms of the effective chemical potentials, v B , by 

= "B + g^B^o + g P Bt3Bbo - (Y B + 9b) A 4 sfn 2 ±6> . (51) 

Then, the zero-temperature limit of Eq. (|32^) yields the baryon grand potential Ob = 
— TlnZff. (We consider only the mean field case here, so AE = 0.) 

The total hadron grand potential 0t o t = ^B + &k = —PV, where P is the total 
pressure. The energy density is given by e = — P + J2b P"B n B + ^K-, where the baryon 
number density, ns, is given by the zero-temperature limit of Eq. (|35| ) and the kaon number 
density is 

n K = ~^ = f(fism 2 9 + 4bsm 2 ±9) with 6 = ]T(Y B + g B )n B /(4/ 2 ) . (52) 

The meson fields are obtained by extremizing Otot; yielding Eq. ( |3"6| ) with T = 0. The 
condensate amplitude, 9, is also found by extremizing Otot- This yields the solutions 6 = 
(no condensate), or, if a condensate exists, the equation 



m s n s B 



fj 2 cos 9 + 2fjb — m 2 K — d\ — d 2 = where 
V 2 di= J2 [{ai + a 2 ){l + Y B q B ) + {ax-a 2 ){q B -Y B ) + Aa 

B^H U H 2 

2/ 2 d 2 sin 2 ±0 = J2 (M* B + gaAa)n B - M A n s Hi - Mxn s H2 

B=H 1 ,H 2 

+ TT^ {Ma ~ Ms)(n ^ " • (53) 
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Here, we have taken g a \ = g a s for simplicity, and the baryon scalar density is 



k F i 



n B = — 



dk k 



{k 2 + M* 2 )l ' 



(54) 



with kpB denoting the baryon Fermi momentum. Equation fl53| ) is equivalent to the re- 
quirement that p, be equal to the energy of the K~ zero-momentum state. 

Finally, we need to satisfy the charge neutrality condition of Eq. (^), which reads 



B 



n K 



0, 



(55) 



and the chemical equilibrium conditions of Eq. (|37|). In the latter, we replace pt^, and 
by nm and /^h 2 ; thus, the first of Eqs. ( |37| ) becomes fijji = fJ-H 2 = Pz° = Mn- Also, 
chemical equilibrium in the reaction n *-* p + K~ requires that the kaon chemical potential 
satisfy pt = p n - p p . 

We shall also need the optical potential for a kaon in nuclear matter for comparison 
with the meson exchange approach. We can use Eq. (^) with equal numbers of neutrons 
and protons, no hyperons, and no condensate (x = 0). Lagrange's equations for an s-wave 
K~, with K- = k-{-x)e- iEt and E = 



p 2 + rn 2 K , can be written 



\V 2 + E 2 



m 2 K ]k 



3nE E KN n s 



4/ 2 



k~U) 



2m K U^ pt k (x) 



(56) 



where n is the density of nuclear matter and n s is the scalar density. For a zero-momentum 
K~ meson, the optical potential reduces to 



U, 



K 



opt 



aK ,yK 
D opt ' v r 



opt 



nK 
D opt 
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2m K f< 



V, 



K 

opt 



3n 



(57) 



Equilibrium nuclear matter density fixes 7£ = —51 MeV, and, for our range of values of 

the sigma term, S* t = -(22-69) MeV. Thus, U* t 100 MeV. This is about half of the 

favored "deep solution" obtained by Friedman, Gal and Batty in their analysis of the 
kaonic atom data, although it is comparable to the value obtained in their t e //P approxi- 
mation. There are, of course, uncertainties in the analysis and also in simply expropriating 
the real part of a complex potential. 
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4.4.2 Meson exchange formalism 

In this approach, we take a Lagrangian for the kaon sector, which contains the usual kinetic 
energy and mass terms, along with the meson interactions, 

C K = d^K+d^K" - {m\ - g aK m K a)K + K- 

+i + 9 P K^\ (K + d, t K- - K~d li K+) . (58) 

Here, b^ 1 denotes the p° field. (The vacuum kaon mass, rriKi is present in the third term 



to render g a x dimensionless.) Schaffner and Mishustin [43] have included an additional 
four-point interaction in their Lagrangian, so that Lagrange's equations yield d^uj^ = as 
required for a particle of spin 1 [^6| . We would argue that at the mean field level the vector 
fields are constants, so the divergence is necessarily zero, and Eq. ( |58| ) is sufficient. We can 
simplify notation by introducing an effective kaon mass defined by 

m* K 2 = m 2 K - g a Km K a . (59) 

Since only the time components of the vector fields survive, it is also useful to define 

X = gwKUJo + g P Kbo. (60) 

In order to determine the kaon partition function at finite temperature, we generalize 



the procedure outlined in Kapusta [13]. First, by studying the invariance of the Lagrangian 
under the transformation — > K^e^ 01 , the conserved current density can be identified. 
The zero component, i.e. the charge density, is 

J = i(K + d K- - K~d K + ) + 2XK + K~ . (61) 
Next, we transform to real fields <pi and fa, 

K± = (fa ± fa)/V2 , (62) 
and determine the conjugate momenta 

tti = d fa - Xfa ; 7r 2 = d fa + Xfa . (63) 

The Hamiltonian density is TLk = ^idofa + ^dofa — £*k, and the partition function of the 
grand canonical ensemble can then be written as the functional integral 



Z K = I [dm][dir 2 ] / [dfa][dfa] 

'periodic 

13 



xexpjy dr J d 3 x (ini-^- + i7T 2 -^- - Hicifa, ^i) + fJ-Jo(fa, 
o 



(64) 
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Here, /x is the chemical potential associated with the conserved charge density, and the fields 
obey periodic boundary conditions in the imaginary time r = it, namely 0$(x, 0) = <f>i(x, f3). 

The Gaussian integral over momenta in Eq. ( ]64| ) is easily performed. Next the fields 
are Fourier decomposed according to 



H =f9co B a + M J £ e t(p - x+ "" T Vi,n( P 

V n,p 

b 2 = f9sma + J yYl 

V n.D 



2,n(p) , 



(65) 



where the first term describes the condensate, so that in the second term 01,0 (p = 0) = 
02,o (p = 0) = 0. The Matsubara frequency oj n = 2imT. The partition function can then 
be written 



Z K = N ( n[#i,«(p)][#2,n(p)]e 5 , where 

ra,p 



S = \pVf 2 9 2 ^ 2 + 2fiX - mf) - § -nC-p), 2 ,-n(-p))D 

n,p 



Wp) 

^2,n(p) 



2 /^+/ + mf -2/iX-/i 



We define the if ^ energies according to 



2(/i + Z)^ 



u# + p z + mf - 2fiX - n 



(66) 



(67) 

and, in most cases, we shall suppress the explicit dependence of oj^ on p. With this 
definition, the determinant of D is 



det D = /3 u4 + (aT - m) ^ + {lo + + nY 



Then 



In Z K = \f3Vf 2 6 2 (iJL 2 + 2/j,X - mf ) - \ ^ln detD . 

n,p 



(68) 



(69) 



Performing the sum over n and neglecting the zero-point contribution, which is not appro- 
priate to a mean field approach, we obtain 



In Z K = yVf 2 9 2 (fi 2 + 2^X - mf ^ 



-V 



d 3 p 
(2^f 



ln(l - e 



+ ln(l - e 



(70) 
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The partition function for the baryon sector takes precisely the form given in Eq. 
and the total partition function is Z tota \ = Z b Zk- Extremization of Z tota \ yields the fields. 
It is useful to define the thermal quantities 



A 



TH 



(2tt) 3 

(ftp 



(p 2 + mf + X 2 ) 1 [f B (u- - //) + f B {u + + (j,)] , (71) 



K J (2vr)3 

where the Bose occupation probability f B (x) = ieP x — 1) _1 . Then the field equations are 



1 ■ -Y,9u,Bn B -gu 1 K[{fe) 2 V + n T K H -XA T K H 
B 

m l b o = 9 P Bh B n B - g P K \{f0) 2 n + n T K H - XA T K H 



B 



m 2 o 



da + ^ 9aB J (2vr)3 V 6 + L J 

+\9a K m K [{fe) 2 + A T K H ] . (72) 



Notice that the condensate contributes directly to the equations of motion (|72|), whereas 
in chiral models the contribution appears in the effective chemical potentials and effective 
masses. The condensate amplitude, 9, is also found by extremization of ^total- This yields 
the solutions 8 = (no condensate), or, if a condensate exists, the equation 

[i 2 + 2fxX -rri$ = \pi-w- (0)] \jm + uj + (0)] = . (73) 

Since fi is positive here, we only have the possibility of a K~ condensate when /i = oj~ (0). 
We illustrate the behavior of the kaon energies with density in Fig. 17. As the chemical 
potential [i increases, a condensate forms when ji becomes equal to u~(Q). At higher 
densities, a condensate is still present so that \x remains equal to lo~(0). Whether the 
tj + (0) energy increases rapidly with density, as in Fig. 17, or remains more nearly constant, 
depends on the strength of the couplings employed. Notice that, utilizing Eq. (|72] ) in Eq. 
(|73|), one obtains a threshold (9 infinitesimal) equation resembling Eq. (|53|) of the chiral 
case; although the weightings of the various baryons are different, and the dU /da term does 
not play a role in the chiral case (see the further discussion below). Above threshold, 9 
enters in different ways in the two models. 

The baryon thermodynamic variables are given by Eq. (|35|). For the kaons, the partition 
function, Eq. fl7C|), gives the pressure, Pk = TV^ 1 InZx; notice that Eq. (|73|) indicates 
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that the condensate gives zero contribution to the pressure. The kaon number density is 
easily found to be 

n K = f 2 e 2 (^ + X) + nl H . (74) 



If we write the baryon energy density in the form of Eq. (p5|), then, after using the equations 
of motion Eqs. (|72| ) and (|7lf), the remaining part of the energy density, which arises from 
the kaons, can be written 

e K = (f0m* K ) 2 + J ^3[ w -(p)/ B ( w --/,)+a; + (p)/ B (a; + +/,)]+X[nr-^r]-(75) 

The total entropy density can be obtained from the standard thermodynamic identity 

s = 0(e + P - [xn K - ^2 hbub) ■ (76) 

B 

The above equations can be applied at zero temperature, in which case the thermal Bose 
occupation probabilities are zero, and the Fermi occupation probabilities become step func- 
tions cut off at the Fermi momentum. 

As before, the neutron star must be charge neutral, i.e. 

^2 qBns - uk - n e - = , (77) 

B 

and in chemical equilibrium, with Eq. (|37]) satisfied and \i = \x n — [i v . 

Finally, as in the chiral model, we can determine the value of the optical potential felt 
by a single kaon at zero momentum. The analogue of Eq. (|5q ) gives 

u o%t = So P t + Vopt ; Sopt = -\g<TK<* ; = -g^K^o ■ (78) 

Since we want to compare the chiral and meson exchange approaches, we will demand that 
they yield the same optical potential in nuclear matter. This should be a reasonable way 
of ensuring that the parameterizations are compatible. Thus, we choose S^ pt and V^ pt to 
be equal to S~p t and 7£ from Eq. (|57|) . For the as yet undetermined kaon-rho meson 
coupling, we take gpx/dpN = 1/3, as suggested by naive quark counting. 



4.4.3 Kaon condensation in non-strange baryonic matter 

Let us begin by studying the critical, or threshold, density at which kaons start to condense. 
The equations obtained in both the chiral and the meson-exchange models can be written 
in the form 

fi 2 + 2/ia - m*K = . (79) 
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For the present, the only baryons we consider are nucleons, in which case the chiral model 
expressions for a and m* K are 



2n p + n r 



2/2 

*2 2 i 

m K = m K + 



2ain s + (2a 2 + 4a 3 )(n,p + n s r 



m s 
2P' 



(80) 



respectively, in nucleons-only matter. In the meson-exchange model, at zero temperature, 
we have 



a 



*2 



(Grn ~ \G p KN )n n + (G KN + \G P KN 
1 dU{a) 



m k + G a KN m K 



9<tN da 



n: 



KN 



where we have used the definitions G 
to isospin symmetric nuclear matter with \i 
the form 



giNgiK/mf, with i 



KN 



m\ = G 



KN 



n 



and 
1 dU(a) 



g aN da 



(81) 

<7, ijj and p. If we specialize 
0, Eqs. © and (|l|) take 



(82) 



respectively. These results give the critical density for condensation in symmetric nuclear 
matter to be u > 6.5, well in excess of the values of u cr n shown in Table 13 for neutrino-free 



matter in beta equilibrium. Schaffner et al. [47] have recently emphasized this and shown 
that, depending on the chosen parameters, condensation may not occur at all in nuclear 
matter. 

The presence of leptons in stellar matter lowers the critical density for condensation by 
a significant amount. In Table 13, we list the critical density for kaon condensation at T = 
for three choices of the constant a^m s , and for matter without and with hyperons. (The 
latter case will be discussed in the next subsection.) Results shown include the neutrino- 
free and trapped neutrino cases. In all cases shown here, baryons are described using the 
GM model, but kaons are described using both the chiral and meson exchange models. 

Consider first the zero temperature case, in which neutrinos are absent. The critical 
densities show a marked reduction as the magnitude of 03771, s increases, since this enhances 
the interactions; but they are not very sensitive to the choice of the model used to describe 
kaonic interactions (i.e., chiral versus the meson exchange model) as long as compatability 
of the kaon optical potentials is required. The mean field model, GM, yields somewhat 



lower values for the critical density than the MRHA models (see Ref. [41]), due to the more 
rapid increase of the scalar fields that enter the interaction terms of Eqs. ( |80|) and (|8l|). 



35 



We remark that if we were to replace the scalar density in these calculations by the 
number density, the critical density ratio, u cr it, for condensation at T = drops by approx- 
imately 1 unit in the density ratio u. (The effect appears to be of smaller magnitude in the 
recent work of Maruyama et al. [38] employing mean field theory.) 

For neutrino- free matter, the effects of finite temperature on the onset of kaon con- 
densation are shown in Fig. 18. (The meson exchange formalism is used for the finite 
temperature results in this subsection.) The results shown are for the GM model and for 
three different values of Y< KN . The proton and kaon concentrations in the top two panels, 
and the chemical potential /x in the third panel, refer to the values at the critical density 
ratio u cr it shown in the bottom panel. Thermal effects give rise to a non-negligible net 
negative kaon concentration, which results in larger proton concentrations than those for 
the zero entropy case. This hinders the onset of condensation. However, the compensat- 
ing changes in the chemical potential \x with temperature (in the range that supports an 
entropy per baryon up to 2) result in only small net changes in u cr n for condensation. Sim- 
ilar results are obtained in the MRHA models. Thus, the effects of condensation remain 
significant, even at finite entropy. 

In Fig. 19, the relative concentrations, the electron chemical potential, and the hadronic 
and leptonic contributions to the entropy are shown for S = 1 for the case of neutrino-free 
matter. The results are for the GM model with a?,m s = —222 MeV. At finite temperatures 
and at densities below that for condensation, the kaon concentration remains smaller than 
that of the other particles; consequently, the changes induced in the total pressure and 
energy density are small. As in the case of zero temperature, the kaon concentration builds 
up rapidly for densities above the condensation density, which results in a softening of the 
equation of state. For an entropy per baryon up to 2, the relative concentrations of the 
various particles essentially retain the zero temperature behavior. Thus, thermal effects 
do not change drastically the softening induced by the condensation. The changes in the 
structure of the star, in particular, the changes in the maximum mass, are therefore at the 
few per cent level, as in the case of stars without condensation. 

We turn now to the effects of trapped neutrinos. The results for this case are contained 
in Table 13 and Figs. 20 and 21. Due to the behavior of the chemical potential (which 
is similar to that shown in Figs. 6 and 10), the critical density for kaon condensation is 
much higher when neutrinos are trapped than in the case of neutrinos having left the star. 
(Compare the critical densities with the central densities in Table 9.) Thus, in the trapped 
case, the hadronic pressures are relatively larger for a substantial range of densities than in 
the neutrino-free case. 

In Table 14, we give the stellar properties for the neutrino-free and neutrino-trapped 
cases. Comparison of the neutrino-free results with Table 8 shows rather small effects, 
except for the largest magnitudes for the parameter a^m s . It is only for these cases that 
the central densities significantly exceed the critical densities (Table 13) and allow a size- 
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able core of condensed kaons to appear. This also occurs in the MRHA models for the 
larger magnitudes of a^m s . Turning to the trapped case, we see that the maximum mass 
is generally larger. This can be contrasted with Tables 8 and 9, where K~ particles were 
absent; and the maximum mass was a little less in the trapped case. This qualitative change 
engendered by kaons is similar to that previously noted for hyperons. 



4.4.4 Kaon condensation in strangeness-rich baryonic matter 

We have affirmed the importance of hyperons in neutron stars, so it is interesting to see 



how they affect the phenomenon of kaon condensation [41|. The necessary formalism has 
been outlined above. Since finite temperature effects are not too large, we focus on the zero 
temperature case. For the chiral model with a^m s = —222 MeV in conjunction with the 
mean field GM description of the baryons, the results are shown in Fig. 22. The critical 
densities are given in Table 13. These critical densities are higher than those in the case in 
which hyperons were absent. The reason is clear. Once a significant number of negatively 
charged hyperons are present (panel (a) of Fig. 20), the electron chemical potential, jj,, 
begins to decrease with density (panel (b)). Since kaon condensation occurs when this 
same chemical potential equals the energy of the zero-momentum K~ state, it is necessary 
to go to higher density where the interactions are able reduce the energy further. For the 
MRHA, in almost all cases, the critical density is beyond the central densities in Table 11, 
so that condensation does not occur. On the other hand, it can take place in the mean field 
model, GM, provided the magnitude of a?,m s is not too small. This follows from the fact 
that the scalar densities increase more rapidly in mean field models than in the MRHA, 



which affects the interaction terms d\ and d 2 in Eqs. (53). Some insight into the role of the 
scalar densities may be gained by examining the threshold condition in the chiral approach 
when only S~ and A hyperons are present: 



2 (2n p + n n - n s -) 2 

M + " V~ m K 



2<nn s p + (2a 2 + 4a 3 )(nj; + < + n|_) 

+ (f(ai + a 2 )+4a 3 )r4] S = 0. (83) 

The first two terms in this equation are smaller than in the nucleons-only case, and this 
has to be compensated by the last term, which requires a higher density. 

The presence of hyperons causes the condensate amplitude to increase rapidly with 
density (panel (d)); so rapidly, in fact, that large changes are induced in the scalar densi- 
ties (panel(c)) and the Dirac effective masses of all the particles (panel(b)). The nucleon 
effective masses (see Eq. (|50[)) are given by 

M* = M - g an a + (2a x + 2a 2 + 4a 3 )m s sin 2 ±6> , 



37 



= M — g an a + (2a 2 + 4a 3 )m s sin 2 ±# , 



(84) 



and go to zero before the central stellar density is reached. This indicates the need to 
consider improvements, in particular, an exact evaluation of the zero-point energy with the 
non-linear Kaplan-Nelson Lagrangian, and this is under investigation. 

Summarizing, we find that while the effect of non-zero temperature upon the onset of 
kaon condensation is small, both the presence of hyperons and neutrino trapping inhibit 
condensation, although the latter is a transient effect. It appears that in the MRHA, at 
least as formulated here, kaons play a much smaller role than suggested by non-relativistic 
treatments [||. Mean field approximations give a larger effect, but we can not yet treat 
them satisfactorily when hyperons are present. 



4.4.5 Sensitivity of kaon condensation to hyperon couplings 

In the calculations above, we assumed that the couplings of the £ and H were equal to 
those of the A hyperon. Here, we relax this assumption and explore the sensitivity to 
unequal couplings of the different hyperons using the meson exchange formalism at zero 
temperature. Of the many possibilities, we pick three for study. These are listed in Table 
15, in terms of the ratio to the nucleon couplings as defined earlier. For the A, we use the 
values discussed previously. For the S, we use two sets of values which gave satisfactory 



fits to the S~ atom data in the work of Mares et al. [48]. This was based on a mean field 
description of nuclear matter using the nucleon couplings of Horowitz and Serot [49], who 
did not include non-linear terms (U(a) = 0). The parameters for this model, termed 'HS81' 
here, are: 

^ = 3.974 fm, ^ = 3.477 fm, and ^ = 2.069 fm. (85) 

Partly for consistency and partly because this model is often used as a baseline in the 
literature, we will adopt these parameters. (Qualitatively similar results are obtained for 
other values of the nucleon couplings, which yield more realistic values of the compression 
modulus.) Finally, we need the couplings of the S. Since there is little information, we take 
the couplings to be equal to those of either the A or the X. Note that case 1 in this table 
is close to the set that we have been using in the previous discussion. 

In Fig. 23, the particle fractions shown in the upper, center and lower panels refer, 
respectively, to hyperon coupling cases 1, 2 and 3 of Table 15. The upper panel is similar 
to results already discussed; note that kaons do not condense up to the maximum density 
displayed, u = 4.5. In discussing the other cases, we first mention the seeming paradox that 
increasing the coupling constants of a hyperon species delays its appearance to a higher 
density. The explanation |3^, [3l]] is that the threshold equation receives contributions from 
the cr, u> and p mesons, the net result being positive due to the uj. Thus, if all the couplings 



38 



are scaled up, the positive contribution becomes larger, and the appearance of the particle 
is delayed to a higher density. With this in mind, consider the center panel of Fig. 23, which 
corresponds to case 2 of Table 15. The £ couplings are larger than in case 1 (upper panel), 
so the £~ no longer appears, thus allowing the chemical potential fi to continue rising with 
density. This allows the H~ to appear at u = 2.2, essentially substituting for the £~ . Of 
course, were we to reduce the E couplings on the grounds that this hyperon contains two 
strange quarks, the would appear at an even lower density. Turning to the lower panel 
of this figure, we recall that this corresponds to case 3 of Table 15, for which both the £ and 
H couplings are increased. Neither of them now appear, and since the chemical potential, 
fj,, continues to increase with density, it becomes favorable for kaons to condense at u = 3.6; 
the fraction Y^-, however, remains rather small. 

Clearly, the lesson to be drawn from this is that the thresholds for the strange particles, 
hyperons and kaons, are sensitive to coupling constants that are poorly known. In matter 
where hyperons are allowed to be present, generally the effects of kaons are small. In fact, 
Schaffner and Mishustin [43 1 find that, with their choice of coupling constants, kaons do not 
condense. On the other hand, should it turn out that the coupling constants of the S and 
3 are larger than the ones adopted here as the standard choice, these hyperons might not 
be present at all, and consequently kaons would play a more important role. Thus, while 
strangeness plays a significant role in determining the constitution and physical properties 
of a neutron star, the detailed behavior cannot be tied down at the present time. 



4.4.6 Metastability of neutron stars with kaon condensates 

Fig. 24 shows the window of metastability in the baryonic mass, which, in the absence of 
mass accretion, is unchanged during the evolution of the star. Here the baryons are nucleons 
described in the mean field GM model without and with kaons, for which we use the chiral 
model with a^m s = —222 MeV. When kaons are present, the range Mb = 2.09— 2.15M Q can 
be supported by the initial EOS of lepton-rich matter, but not by the later EOS of lepton- 
poor matter (lines ending in dots). This range of metastability corresponds to gravitational 
masses of Mq = 1.81 — 1.91-Mq. In the absence of kaons, metastability does not occur, 
since the maximum mass decreases when the neutrinos leave (lines ending in stars). Note 
that the qualitative features here are similar to the case of matter with strangeness-bearing 
hyperons, see Fig. 15. 

Fig. 25 shows the corresponding evolution of the maximum gravitational mass in mat- 
ter with and without kaons in the deleptonization stage. As we have previously remarked, 
the initial state with Yl b = 0.4 corresponds to an electron-neutrino fraction of Y Ue ~ 0.08 
and, of course, in the final state Y Ve = 0. As in Fig. 16, the decrease (increase) in the 
maximum mass when kaons are present (absent) is most pronounced shortly before the 
neutrino fraction drops to zero. 
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4.5 Matter with quarks 

We now examine another type of softening of high density matter by allowing for a hadron 
to quark phase transition in the interior of the star . Glendenning |5(| has suggested that 
a mixed phase of baryon and quark matter exists over a wide range of densities in the case 
of neutrino- free matter. Our primary interest here is in protoneutron stars and the effects of 
trapped neutrinos ||], which, conceivably, could lead to observable consequences. We focus 
on zero temperature, since, as we have seen, changes in the maximum mass due to neutrino 
trapping are larger than those due to finite temperature. The influence of complicated 



finite size structures due to Coulomb and surface effects [51] does not qualitatively affect 
our conclusions and will be taken up elsewhere. 

We shall follow fairly closely the treatment of Glendenning [jSCfl . Thus, for the pure 
phase in which the strongly interacting particles are baryons, we employ the mean field GM 
model, for which the formalism has been discussed in earlier sections. For the pure quark 
phase (in a uniform background of leptons) , we use the bag model for which the pressure is 



kF f ,Q, , O k 



1 V d 3 k k 2 1 ^ "[' d 3 k 

Q -B + - ^ g f J — +-^- 

o 1 1 o 

The first term accounts for the cavity pressure, and the second and third terms give the 

Fermi degeneracy pressures of quarks and leptons, respectively. The constant B has a 

simple interpretation as the thermodynamic potential of the vacuum, and will be regarded 

as a phenomenological parameter in the range (100 - 250) MeV fuT 3 . The lower limit here 

is dictated by the requirement that, at low density, hadronic matter is the preferred phase. 

For B much larger than the upper limit, a transition to matter with quarks never occurs. 

The degeneracy factor for quarks is gf = 2 x 3, accounting for the spin and color degrees 

of freedom. The chemical potential of free quarks in the cavity is fif = ^Jkpj + m 2 , where 
kpf is the Fermi momentum of quarks of flavor /. For numerical calculations, we take the 
u and d quarks as massless, and m s = 150 MeV. The baryon density and the energy density 
are 

l \ -> K Ff 

n Q=3 L "/> U f = ^2 

f=u,d,s 

£q = -Pq + n fVf + nim . (87) 
/ i 

The relevant weak decay processes in the pure quark phase are similar to Eq. (|l|), but with 
Bi replaced by qf, where / runs over the quark flavors u , d, and s. In neutrino- free matter, 



3 f _^, 1 J (2vr)3 (fc2 +m 2)l/2 3V 7 (2VT)3 (A;2 +m 2)l/2 ' 
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charge neutrality and chemical equilibrium under the weak processes imply 

£ if n f + Z Q fne = 

fJ-d = flu + M = Ms • (89) 

When neutrinos are trapped, the new chemical equilibrium relation is obtained by the 
replacement \xg_ — » /i£ — /i^ in Eq. p9| ) . 

In the mixed phase of hadrons and quarks, it is necessary to satisfy Gibbs' phase rules: 

Ph = Pq and /J, n = fi u + 2(j, d . (90) 



Further, following Glendenning |50J, we require global, but not local, charge neutrality of 
bulk matter, for both separately conserved charges: baryon number and electric charge. 
Denoting by / the fraction of volume occupied by the hadronic phase, we have 

/5>B71B + (1-/) E <7/?V+E^ = (91) 

B f=u,d,s £=e,/x 

n = /5>B + (1- f)n Q , (92) 

where qs is the electric charge of each hadron. An important consequence of global neu- 
trality is that baryonic and quark matter coexist for a much larger range of pressures than 
for the case of local charge neutrality [50]. The total energy density is e = fen + (1 — f)sg- 



4.5.1 Metastability of neutron stars with quarks 

Fig. 26 shows a comparison of the compositions of neutrino-free matter (top panel) and 
neutrino-trapped matter (bottom panel). In the case of neutrino-free matter, quarks make 
their appearance at around Auq for B = 200 MeV fm -3 . After this, the neutral and 
negative particle abundances begin to fall, since quarks furnish both negative charge and 
baryon number. The bottom panel of Fig. 26 shows the influence of trapped neutrinos 
(with Ylb = 0.4) on the relative fractions. The primary role of trapped neutrinos is to 
increase the proton and electron abundances, which strongly influences the threshold for 
the appearance of hyperons. The A and the S's now appear at densities higher than those 
found in the absence of neutrinos. In addition, the transition to a mixed phase with quarks 
is delayed to about 10no- Qualitatively similar trends are observed for other values of the 
bag constant B. 

In Fig. 27, we show the phase boundaries as a function of the bag pressure B. The 
onset of the phase transition is at density n\ = u%no, and a pure quark phase begins at 
density ri2 = U2Uq. Also shown are the central densities n c = u c n$ of the maximum mass 
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stars. Trapping shifts the onset of the phase transition to higher baryon densities and also 
reduces the extent of the mixed phase in comparison to the case of neutrino-free matter. 
The existence of a mixed phase inside the star depends on whether or not hyperons are 
present. In the absence of hyperons (top panel), a mixed phase is present for the entire 
range of bag constants. When hyperons are present (bottom panel), the mixed phase is 
present only when B is sufficiently low (B < 165 MeV fm~ 3 ) and occurs over a smaller 
range in density than that found in the absence of hyperons. The abrupt change in the 
onset of the transition around B = 140 MeV fm~ 3 is caused by the appearance of hyperons 
prior to that of quarks. 

The dashed lines in these figures correspond to the case in which neutrinos have left 
the star. Whether or not hyperons are present, the mixed phase is now present over a wide 
range of density inside the star. Note also that, since the central density of the star u c < U2 
for all cases considered, the presence of a pure quark phase is precluded. Finally, note that 
quarks are more likely to appear the fewer the neutrinos remaining in the star. 

Table 16 shows the maximum masses of stars as a function of the composition of the 
matter. With only nucleons and leptons (last row), neutrino trapping generally reduces 
the maximum mass from the case of neutrino- free matter. This is caused by the smaller 
pressure support of lepton-rich matter, in which the gain in the negative symmetry pres- 
sure exceeds the increase in leptonic pressure. However, the introduction of quarks, which 
soften the EOS, causes the maximum mass for the trapped case to be larger than that for 
neutrino-free matter. This reversal in behavior is due to the fact that the first appear- 
ance of quarks occurs at a higher density when neutrinos are trapped. When hyperons are 
present, the maximum mass remains larger for the trapped case. This is also due to the 
appearance of hyperons at a higher density when neutrinos are trapped. It is an example of 
the general result that when matter contains non-leptonic negative charges, the maximum 
mass of the neutrino-trapped star is larger than that of the neutrino-free star. This result 
has important ramifications for the evolution of proto-neutron stars and for the formation 
of black holes, as we discuss in Sec. 6. 

4.6 Global energetics 

While the softening of the EOS due to the presence of negatively-charged particles has 
a large effect upon the maximum mass, it has surprisingly little effect upon the binding 
energy versus mass relationship for neutron stars. The binding energy is the difference 
between baryonic and gravitational masses of the final neutron star configuration. It is an 
important observational parameter, because at least 99% of it appears as radiated neutrino 
energy. In Fig. 28, we display the binding energy as a function of the baryonic mass Mb 
for stars with and without the presence of strangeness-bearing components, such as kaons 
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or hyperons. The various curves refer to different EOSs and terminate at the maximum 
mass supported by their respective EOS. 



A few striking results are evident from this figure. 

1. The largest binding energy occurs for the EOS that supports the largest maximum 
mass. 

2. For each EOS, the binding energy displays a nearly quadratic behavior up to the 
maximum mass. 

3. There exists a rather narrow band of possible binding energies for a given mass, 
implying the following universal relationship for the binding energy as a function of 
mass: 

Mb\ 



B.E. = (M B - M G y * (0.065 ± 0.01) ( ^ ) M Q , (93) 



where the numerical coefficient represents an update of the value quoted earlier by 
Lattimer and Yahil . The universality is not altered by the presence of significant 
softening in the high density EOS due to the appearance of quarks, kaons, or hyperons. 

4. Only near the terminations at the maximum masses do the binding energies slightly 
deviate from the lower envelope of the curves. This effect is slightly more pronounced 
for softer equations of state. 

We have found that the lower envelope of the binding energy-mass relation is equivalent 
to that found for the stiffest plausible equation of state, namely one that is limited by 
causality. Denoting by n t a transition density above which the EOS is assumed to be 
causal, the EOS above n t is given by |53]|: 



Pt ~ e t + (Pt + e t ) - 
n t 



n x 2 



1 

2 

et + P-Pt, (94) 



where the quantities Pt and et are the pressure and energy density at nt and thus depend 
both upon it and upon the equation of state employed. However, if nt is in the range 
no — 2no, both Pt and et are somewhat insensitive to the EOS, and the binding energy is 
relatively insensitive to these quantities. Fig. 28 includes the binding energy as a function 
of baryon mass for an EOS which is causal above nt = 0.3 fm~ 3 and matched to the GM 
equation of state below nt- Varying the values of nt and the EOS below m merely alters the 
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termination point (maximum mass) of this curve, without otherwise noticeably changing 
it. 

We note that an analytic representation of the binding energy-mass relation can be 



determined using the techniques developed by Nauenberg and Chapline [54|, who assumed 



that the pressure and energy density are constant within the neutron star as an alternative 



to an explicit integration of the Tolman-Oppenheimer-Volkov equation [14]. Defining the 
parametric variable x by 

the mass is 



/ 3c 6 (l-£) , , 6cosx , v 

y 327rG" 5 (e f - Pt) 9 cos x ~ Q 



where Q = 2sin 3 x/(x — sin x cos x)- The binding energy can then be written as 

B.E. = M \ * UtmBc2 v/I ^ ? - ll . (97) 

Utilizing Q = 3(1 — 3x 2 /10 + • • •) and ^ = x 2 /10 + • • •, as appropriate for low mass stars, 
one recovers the Newtonian result that B.E. = (3/5)GM 2 /Rc 2 + • • •, which displays the 
quadratic dependence on the mass of the star. 

Since nearly all of the binding energy is released in the form of neutrinos, it appears 
that an accurate measurement of the total radiated neutrino energy will lead to a good 
estimate of the remnant mass. However, as the results of Fig. 28 and Eq. ( |9"3| ) show, it will 
not be possible to distinguish the various equations of state from the total binding energy 
alone. 



5 Evolution timescales 

We have seen that the structure of a neutron star is strongly influenced by the presence of 
trapped neutrinos and, to a lesser extent, by the non-zero entropy/baryon. The evolution 
of the star will be governed by the timescale for release of the trapped neutrinos, referred 
to as deleptonization, and the timescale for thermal cooling, which reduces the entropy to 
a small value. We discussed these timescales in Sec. 1 on the basis of detailed numerical 
calculations of the evolution dynamics ||]. Since the timescales are of some significance, 
we will show in this section how these can be estimated from analytical considerations and 
illustrate the dependence on the equation of state and the opacities. 
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The equations that describe the physical state and evolution of the nascent neutron 
star are simply those of standard stellar structure theory, modified for the effects of general 
relativity and augmented to include lepton and neutrino transport. (Photon transport is 
completely suppressed at the high densities of the core.) All the known neutrino species 
and their antiparticles carry energy. Most of the time, in most of the star, the neutrinos 
are in thermal equilibrium with the matter and have Fermi distributions. This is true 
because neutrino processes, such as the nucleon Urea, pair production, and inverse nucleon 
brehmsstrahlung, are sufficiently rapid. The calculations of Maxwell [E^], for example, 
illustrate that the timescale for v e {yS) equilibration will be less than 1 s, at nuclear densities, 
for temperatures T > 2(3) MeV. This is an overestimate for the equilibration time for u e , 
since Maxwell did not consider the direct Urea process. If that process is included, the 
timescale decreases by an order of magnitude. 

It is a good approximation to lump together and v T transport together into "yu" 
transport, and thereby assume them to have zero chemical potential and equal opacities. 
Further simplification can be made if one assumes the electron neutrino and antineutrino 
opacities are also equal. Then we need only define one e-type chemical potential, denoted by 
[i v . This is exact deep in the opaque interior, but breaks down in the transparent regime 
above the neutrinosphere. It can be joined smoothly to a free-streaming approximation 
here, however. The relevant equations, in spherical symmetry, have been given in Ref. 
to which we refer the reader for more elaborate discussion. They are 

dP G(M + Airr 3 P)(p + P/c 2 ) 



(98) 

4vrr 2 p (99) 

(100) 



dr r(r - 2GM/c 2 

dM 2 
dr 

dN 47rr 2 n 

~dr ~ y/l - 2GM/rc 2= 
dY v M^r 2 F u e 



dN +3 » < 101 > 
dY 

= -S v (102) 

dU d(l/n) 2^dL u e 2 ^ , , 

z = -P y - e~ 2<i> " . 103 

dr dr ON y ' 



Here Eq. ( |98[ ) is the general relativistic equation for hydrostatic equilibrium, in which M(r) 
is the enclosed gravitational mass. The enclosed baryon mass, N(r), obeys Eq. ( [LOOj ). Eqs. 
(p0l| ) and ( gOp give the rate of change of the electron neutrino and electron concentrations, 
with F v the number flux of electron neutrinos and S v the electron neutrino source term. 
Finally, Eq. ( [L03D gives the rate of change of U, the internal energy per baryon, where L v 
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is the total neutrino luminosity (including all species). The term = \/—goo relates time 
at infinity r with the coordinate time t, and one can show [|J that dcj)/dP = — (P + /jc 2 ) -1 . 

In the diffusion approximation, fluxes are driven by density gradients. In our context, 
this translates into expressions of the form 

K _ .J-^^fA dEp , (10 4) 

L. - "f^Ef (10=) 

where the sum is over neutrino species. The \ v and A^'s are mean free paths for number 
and energy transport, respectively, and are functions of neutrino energy E v . Also, n u (E v ) 
is the number of electron neutrinos with energy E v , and €i(E u ) is the energy density of 
species i = e,fj,. The general relativistic corrections have been dropped for clarity, although 
they are straightforward to incorporate. 



We can combine Eqs. ( 101 ) and ( 102 ) to obtain the rate of change of the total lepton 
number, and Eq. (|103| ) and the first law of thermodynamics to obtain the rate of change of 
the entropy: 

i=n,p,e,v 

To proceed, we have to understand the energy dependence of the opacities. There are 
two main sources of opacity: 

1. Neutrino- nucleon absorption. This affects v e and v e only, except at very high densities 
if muons are present, which is not possible until relatively late in the cooling. The 
absorption mean free path, assuming nondegenerate nucleons, is 

EL ^ 2 



Kbs = Kbs\ 1 ^) cm, (108) 

where \° bs is the fiducial absorption mean free path at the fiducial v e energy E vo ~ 260 
MeV, the typical v e chemical potential at the beginning of deleptonization. From 



Ref. [56 1, we find 

n n a (l + 3g A ) \ E 



K bs = i^rA^\ , (109) 
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where n n is the neutron number density, a Q = 1.76 x 10~ 44 cm 2 , and qa — 1-257. 
Using n n = (8/3) no, appropriate for Y e = 1/3 and a total baryon density n = 4no 
at the beginning of deleptonization, we find X° ahs ~ 0.36 cm. Since nucleons will, in 
fact, be at least partially degenerate, and because of fermi liquid effects, the true 
absorption mean free path will be about 3-10 times larger than this value. 

2. Neutrino-nucleon scattering. This elastic scattering affects all z/-types. For nonde- 



generate nucleons [56|, 



Ks 




(110) 
(111) 



for v e and v^, respectively. For our reference density, 4no, we obtain X° es ~ 1.37 
cm for E vo = 260 MeV, and A° s ~ 1.75 cm for E Vfl0 = 230 MeV; the latter is the 
appropriate value for the mean Vn energy, with s = 2, at the beginning of the cooling 
era. Corrections to the scattering mean free paths for degeneracy and interactions 
should be similar to those for absorption. 

Thus, during deleptonization, \° abs < X°, and v e absorption dominates both energy and 
lepton number transport. However, during thermal cooling, energy transport is effected 
mostly by jx— and r— neutrinos, since f e 's are more tightly coupled to the matter. 
The opacities imply three kinds of fluxes: 



1. a number flux F v of f e 's, dominated by absorption. Utilizing Eq. ( 108 ), we have 



cX abs dn u (E u 
3 dr 



[ dEv = cX °abs E lo g/^ _ J^v 



67r 2 (7ic) 3 dr dr 



;ii2) 



2. an energy flux L e v of z^ e 's, also dominated by absorption: 



2 c\ abs de u {E u ) 

Airr — dE u 

6 Or 



inr a— 
or 




(113) 



3. an energy flux of f^'s and zVs, dominated by scattering: 



3 Or 



4irr z 



4irr 2 b-^- ( - — 
or \ 6 



cKs (Eu.of d /vr 2 T 2 



37r 2 (/ic) 3 dr \ 6 



(114) 
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where we used the fact that v^s and i/ r 's have zero chemical potential. 
5.1 Deleptonization Era 

We can now appreciate the separate stages of deleptonization and cooling. Deleptonization 
is dominated by number transport, and the controlling equation is 

dYi, 8Yl dY v a d f 2 d^v 



n —±L = n ll±L ZUL = r i~rv /hjqx 

dt dY u dt r 2 dr\ dr J ' v ' 

During deleptonization, the electron neutrinos are degenerate, and to lowest order nY v ~ 
(/i u /fic) 3 /Git 2 . The neutrinos and electrons are in beta equilibrium, and it can be shown 
that dY L /dY u = (dY L /dY u ) (E U0 /fi u ), where {6Y L /8Y V ) ~ 3 for Y uo ~ 0.06. Thus, 

3E "° {dYj ^~dr = V 2 d~r V ^rT) ■ (116) 
We now seek separable solutions of the form //„ = E vo cj)(t)ip(r). We obtain 

3 f^A] ^k d A = ( r 2^t\ = _ a (117) 

\ dY v J a dt r 2 ^ 2 dr \ dr J 

where a is a separation constant. One sees that 

* =i - t/T " Td ^M\- (u8) 

where is the diffusion time. Note that the decay of the neutrino chemical potential 
is approximately linear with time. This result is borne out in more detailed numerical 
calculations. The radial equation is 

,2,2 Id/ 2< 9 V , 



x z ox V ox ' 



which is just that of the Lane-Emden polytrope (57j] of index 2. In the above, R is the 
radius of the star. The solution of this equation satisfies aR 2 ~ 19; and, therefore, 

9R 2 ( R \ 2 . . 

r d ~ ~ 44.3 — - — s. 120 

d I9c\° abs VlOkm/ V ; 

Recalling that degeneracy and fermi liquid corrections will increase the mean free path by a 
factor of 3-10, a deleptonization time of 5-15 s is indicated. This is the correct magnitude 
for the deleptonization time, and shows clearly how it depends upon the equation of state 
through the radius of the star, R, and upon the opacity. 
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Because of the positive temperature gradient and the chemical potential gradient, the 
deleptonization is accompanied by heating in the core. The entropy/baryon rises to the 
value of about 2 before it decreases during the cooling era. The onset of cooling does not 
begin until deleptonization is complete. Rewriting the last term of Eq. (|107l) as 



^ dYi ,dY e dY L 

n,p,e,v 0,1 al al 

we can combine Eqs. fllOlf ) and ( |10^ ) to find 

nT ds = a + bd_ ( 2 9(vr 2 r 2 /6) \ + a fdj^A 2 ^ 
dt r 2 dr \ dr J \ dr ) 

The terms proportional to a are due to electron neutrinos, and the term proportional to b 
is due to the other neutrinos. There is heating or cooling depending on the direction of the 
temperature gradient, but the chemical potential gradients always lead to heating. When 
fi u S> T, the [dix v jdr) 2 term dominates, and we have heating. When \x v — 0, and the 
temperature decreases with radius, cooling occurs. 

5.2 Thermal Cooling Era 

We now turn to the thermal cooling of the protoneutron star, which continues beyond the 
deleptonization era. While the initial entropy per baryon s in the star's interior is about 1, 
after the deleptonization heating is finished the entropy reaches the value of about 2. The 
entropy is dominated by baryons for temperatures less than about 100 MeV. Thus, we may 
write [27 1 

1 771* / Tin \ 

s^2a ed T ; a u = — — — MeV" 1 , (123) 

15 m \n J 

where m* is the effective nucleon mass. Note that the maximum value of the central 
temperature is 

/ m \ / n \ 2 / 3 



MeV. (124) 
We henceforth neglect the density dependence of m* and use m* ~ 0.5m. Notice that the 



estimate of T max in Eq. (124) agrees quite well with our previously tabulated results. 

The cooling is dominated by the ii- and r-neutrinos, since b > a. With these simplifi- 
cations, Eq. (|122j ) becomes 



6a £d n dT 2 1 d ( (125) 



7r 2 (a + b) dt r 2 dr \ dr 
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Once again, we may separate the resulting equation in terms of the time and radial depen- 
dence of the temperature. Writing T = T max ip(r)ip(t), we find 

Ylai^n 1 dcf> 1 d ( 2 <9'0 



r-k- =-«> (126) 



7r 2 (a + b)<fi dt ip 2 r 2 dr \ dr 

where a is a separation constant. The solution of the radial equation is that of an n = 1 
Lane-Emden polytrope, for which the eigenvalue is a = ir 2 /R 2 . The temporal equation has 
the solution 

= exp(-(i-r d )/r c ) (127) 
for t > Td, where 

_ l2a£dnR 2 3s max (E uo /T max )R 2 noR 1 ! 

Tc " vr 4 (a + b) ~ Y uo ^c(X° abs /2 + \%(E^ o/ 'E vo ) 2 ) ' 1 ) 



Note that 



r d W vo ^ T max 1/2 + (\°JKbs)( E ^o/E Uo 



1.7. (129) 



This result is independent of our assumptions regarding s max or T max . This is also of the 
right magnitude to match the numerical calculations, which indicate that r c = (1 — 2)r^. 
Thus, in spite of the fact that the mean free paths that dominate cooling are larger than 
those that dominate deleptonization, the large ratio of the matter's heat capacity to that 
of the neutrinos forces the cooling time to be longer than the deleptonization time. 

At late times, however, the decay of the central entropy or temperature is roughly 
linear with time Q. This feature can be seen to be a result of the increasing degeneracy of 
the star. For degenerate nucleons, the mean free paths have an E~^ dependence. Ignoring 
the absorption contributions, this energy dependence leads to a linear time decay, with a 
time constant 

/ _ s max R 2 I E uo \ (1301 
c 76(ln2)Y, oC A' s \E. ' 1 ] 



where the fiducial scattering mean free path for degenerate nucleons is [56] 




K = Ks-r—^r IP- U 1-6A- • (131) 



For the values we have been assuming, one finds that t' c ~ 11 s. Of course, since this 
result is applicable to the later stages of cooling, the use of a smaller value of E v may 
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be appropriate, which will lead to an increase in this timescale estimate. In any case, one 
expects the experimental behavior of Eq. ( |129| ) to alter to a linear decay after an e-folding 
time. 

It is finally interesting to note the overall sensitivity to the equation of state. For a 

— 1/3 

fixed mass star, employing R oc n c , where n c is the central density, we determine that 
Td oc nj^ and r c cx n° so that for nucleons-only matter, the central density affects the 
timescales rather weakly. However, larger changes in the opacities and radii are trigerred 
by the appearance of negatively charged strongly interacting particles in any form [58], 
suggesting that their onset may significantly affect these timescales. 



6 Implications 

Our main thrust in this work has been to elucidate how the structure of a proto-neutron 
star depends on its composition, which is chiefly determined through the nature of the 
strong interactions at high baryon density. During its early evolution, a neutron star with 
an entropy per baryon of order unity contains neutrinos that are trapped in matter on 
dynamical timescales. After a time of a few tens of seconds, the star achieves its cold, 
catalyzed structure with essentially zero temperature and no trapped neutrinos. The in- 
fluence of finite temperature on the star's structure is dominated by the behavior of the 
baryonic thermal pressures, which are governed by the behavior of the baryonic effective 
masses. Baryonic thermal pressures are proportional to their Landau effective masses, so 
that nuclear models that lead to extremely small effective masses at high density, such as 
Skyrme-type interactions, will generally show substantially larger effects at finite tempera- 
ture than other models. We have shown, however, that the gross reduction of the effective 
mass in the Skyrme case leads to acausal sound speeds in dense matter and should be 
discounted. Thus, finite temperature effects upon the maximum neutron star mass are 
naturally limited. 

In general, however, changes in the maximum mass due to neutrino trapping are larger 
than those due to finite temperatures. These changes depend sensitively on the compo- 
sition of matter, in particular, on the question of whether or not a new component that 
substantially softens matter can appear in the cold, catalyzed star at high density. The new 
components that have been discussed to date include hyperons, a pion or kaon condensate, 
and a transition to quark matter. All these components involve negatively charged non- 
leptonic matter; hence, they appear at lower density in the cold, catalyzed star than in the 
hot, neutrino-trapped star. Since thermal pressure is always positive, the cold, catalyzed 
star always has a higher density than a hot, neutrino-trapped star. Consequently, a cold, 
catalyzed star contains the softening component in a larger proportion of the star's mass 
than the hot, neutrino-trapped star, and this leads to a smaller maximum mass. 
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This behavior is opposite to that found for equations of state containing only nucleons 
and leptons and no additional softening component. In this case, neutrino trapping gen- 
erally reduces the maximum mass from the value found in neutrino free matter; although 
neutrino-trapped matter contains more leptons and more leptonic pressure, it also contains 
more protons and, therefore, less baryonic symmetry pressure. While finite entropy pro- 
vides additional pressure support, the amount of the increase over the zero temperature 
case is generally small, especially if realistic forces are employed. 

It must be emphasized that the maximum mass of the cold catalyzed star still remains 
uncertain, due to the uncertainty in strong interactions at high density. At present, all 
nuclear models can only be effectively constrained at nuclear density and by the condition 
of causality at high density. The resulting uncertainty is evident from the range of possible 
maximum masses predicted by the different models considered in this work. Despite this 
uncertainty, our findings concerning the effects of finite entropy and neutrino-trapping offer 
intriguing possibilities for distinguishing between the different physical states of matter. 
These possibilities include both black hole formation in supernovae and the signature of 
neutrinos to be expected from supernovae, as we now discuss. 



6.1 Black hole formation 

The gravitational collapse of the core of a massive star produces a lepton-rich, neutrino 
trapped, proto-neutron star and an expanding shock wave. Energy losses from dissociation 
of heavy nuclei and neutrinos weaken the shock, preventing a "prompt" explosion. Within 
a few milliseconds, the shock wave stagnates into an accretion shock at a distance of 100- 
200 km from the proto-neutron star. Gandhi and Burrows |5!| have demonstrated that 
the neutrino luminosity of the stellar remnant is able to quasi-statically support the shock 
against the ram pressure of infalling matter, which accretes onto the neutron star. Re- 
cent successful models of gravitational collapse supernovae [50, 61] invoke delayed neutrino 



heating, augmented by convective motions, to power the supernova explosion. 

The explosion appears to occur within \ to a few seconds after the core bounce; and, 
once expansion occurs, accretion onto the proto-neutron star is diminished. Therefore, the 
star can be expected to accumulate nearly all of its baryon number within a few seconds of 
core bounce. 

As we have seen, a nucleons-only star has a maximum mass that grows as neutrinos leak 
out of it. It thus appears unlikely that such a star could form a black hole during the longer- 
term deleptonization era. In this case, a black hole could only be produced immediately 
after bounce or during the short-term accretion stage, when the shock is stagnant. 

However, this is not the case for matter with non-leptonic, negatively charged softening 
components. For matter containing hyperons, kaons, or quarks, the maximum mass of the 
neutrino-trapped star is larger than that of the neutrino-free star. Should the maximum 
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baryon mass of the cold neutrino-free star be close to 1.5Mq, black hole formation could 
occur as the neutrinos diffuse out of the protoneutron star, i.e., during the first 10 seconds 
following bounce. Burrows || has demonstrated that black hole formation should be ac- 
companied by a dramatic cessation of the neutrino signal, since the event horizon invariably 
forms outside the neutrinosphere. Such behavior would be relatively easy to observe from 
a galactic supernova and would suggest that the equation of state produced a metastable 
protoneutron star. 

In order to highlight some of the observable consequences for different compositions 
of high density matter, the gravitational and baryonic masses, and their differences, both 
for cold, catalyzed matter and for hot, neutrino-rich matter, are shown in Fig. 29. Two 
generic compositional cases are displayed: nucleons-only matter (npe~ and matter with 
hyperons (npHe~ fi~). Matter with condensed kaons and with a quark-hadron transition 
has a similar behavior, and is discussed in more detail in Refs. [||, [|. 

In all cases, the gravitational mass for the neutrino-rich cases are greater than for the 
untrapped case for the same baryon mass, reflecting the binding energy released during 
this stage of neutron star formation. This is usually less than half the total binding energy 
released, the remainder being emitted during the prior stage. The prior stage is the period 
between core bounce and the production of the approximately adiabatic (s ~ 2), lepton-rich 
star (shown by the upper solid curves in Fig. 29), and lasts about 1-3 seconds. The binding 
energy emitted during the prior stage is difficult to show in Fig. 29, because during this 
period the neutron star is rapidly accreting mass (i.e., Mb is increasing). The accretion 
should drop substantially, and the value of Mb should approach a limiting value after the 
first few seconds. We note that the binding energy emitted in the deleptonization and 
cooling stage appears to be insensitive to the composition of dense matter, just as the total 
binding energy was found to be a universal function of mass. Given its increase with stellar 
mass, this means that the total energy released is apparently not a good discriminant of 
composition. 

However, Fig. 29 clearly shows the consequence of different compositions for black hole 
formation. Black hole formation can be observed as an abrupt cessation of neutrino signal, 
since the event horizon forms outside of the star's neutrinosphere. For the nucleons-only 
case, black hole formation is unlikely to occur during the deleptonization and cooling stage, 
the one marked by the transition between the upper solid and the lower dashed curves, since 
Mb is approximately constant (or only slightly increasing) . If a black hole were to form from 
a star with this composition, it is much more likely to form during the post-bounce accretion 
stage. This is not true for the other compositional cases. Here, the neutrino-trapped matter 
is always capable of supporting more mass than the cold, catalyzed matter. If the baryon 
mass of the proto-neutron star is near the maximum mass of the cold, catalyzed neutron 
star, as it is, for example, for a 1.5 M star in both the npHe~ \i~ and the npHQe" jjT 
cases, then there exists a range of 0.1 — 0.2 M Q above this mass in which a hot neutrino- 
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trapped star can be stabilized. Thus, with this composition, a black hole could form during 
either stage of neutron star formation; although, given the relatively small value expected 
for the mass of the imploding pre-bounce core, it seems more likely that a black hole would 
form during the later stage. 

A consequence of the potentially long delay of 10-15 seconds between core bounce and 
black hole formation is that black hole formation and total binding energy release are not 
necessarily correlated. The softening in the equation of state marked by the appearance 
of negatively charged hadrons is accompanied by relatively little further binding energy 
release. Thus, the large energy release inferred by the neutrino detections from SN 1987A 
did not imply that the equation of state could not yet change due to the appearance of any 
exotic matter. 

Black hole formation, by cutting off the neutrino luminosity from the protoneutron star, 
would short-circuit the supernova mechanism if an explosion had not already occurred. This 
appears to be more likely for the nucleons-only stars or for stars with relatively large initial 
core masses, i.e., very massive stars with M > 20 — 30 M Q ||62| , pSfl . Note also that these 
scenarios have different implications for nucleosynthesis, since prompt black hole formation 
and a successful supernova explosion, in which newly synthesized nuclei are ejected, may 
be incompatible. 

Based on the calculations of Prakash and Lattimer, Brown and Bethe [JTTJ] proposed 
that a window exists in which neutron stars collapse to black holes during deleptonization 
if hyperons, kaons or quarks are present in neutron stars. This proposal was meant to 
explain the apparent non-existence of a neutron star in the remnant of SN1987A, although 
a neutron star may have temporarily existed some 10-15 s, during which neutrino emission 
was observed. The window naturally exists if negatively charged hadronic matter appears 
after deleptonization. The particular case of SN1987A will be discussed below. 

6.2 Neutrino signals from supernovae 

The composition of a neutron star will also influence the details of the star's neutrino 
emission. We focus on two possible diagnostics - the total radiated energy and the relative 
numbers of emitted neutrinos of different types. 

In Fig. 30, we display the electron concentration, Y e , for various dense matter compo- 
sitions. The reference straight line is the total electron-lepton concentration of the initial 
proto-neutron star, which we have assumed to be approximately Yl £ = 0.4. The difference 
in any of the other curves from the reference line shows, as a function of density, the total 
net electron-neutrino concentration (specifically, the difference of the v e and v e concentra- 
tions) that eventually leaks out of dense matter. The integral of the difference over the 
stellar density profile gives the net number difference for the entire star. This difference is 
always positive, although the actual value is sensitive to composition. 
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A rough measure of the relative fluxes of escaping neutrinos can be found as follows. 
Although, in the stellar core, there are essentially only i/ e 's because n e /T S> 1, fj, e 3> 1 
and fj,a = \i T ~ 0, the emerging neutrino signature consists of nearly equal proportions of 
all six types of neutrinos. This may be understood as follows. Diffusion degrades the high 
energies (fj, Ue ~ 200 MeV) of the core f e 's; therefore, the emerging neutrinos have average 
energies in the range 10-20 MeV. Pair production in the hot matter in the outer mantle 
of the proto- neutron star generates several pairs of all three neutrino flavors per core v e . 
Using 10 MeV for the emergent energy, we find about 3 pairs emerge per core v e emitted. 
This shows that the total number of escaping neutrinos has a slight excess of f e 's; for this 
example, the ratio of v e : v x , where x refers to any of the other neutrino species, is about 
4:3. 

The actual situation is more complicated because, in general, different neutrino species 
are emitted with different energies. Also, the total flux of neutrinos should be sensitive to 
the available binding energy change. Nevertheless, the basic trends are clear: First, the 
smaller the value of Y e in the cold, catalyzed neutron star, the larger the excess of z/ e 's 
that will be emitted. Second, the larger the binding energy change during deleptonization, 
the smaller the relative excess of v e 's will be. This is because the energies of the escaping 
neutrinos are rather insensitive to the structure of the proto-neutron star, including details 
of the equation of state of high density matter. The energies of the escaping neutrinos are 
determined instead by the properties of the outer mantle of the star. Thus, higher binding 
energy release translates directly to larger total numbers of escaping neutrinos. 

However, referring to Fig. 29, one sees that the change in binding energies during delep- 
tonization and cooling is relatively insensitive to the dense matter composition, although 
it does increase with the final neutron star mass. Therefore, the net excess of zVs from 
a proto-neutron star during deleptonization and cooling seems to be a probe of the final 
value of Y e in dense matter, a quantity that is quite sensitive to dense matter composition. 
It remains to be seen if differences in the final value of Y e are large enough to be observable 
in the neutrino signal from a galactic supernova. 



6.3 Supernova SN1987A 

The case of SN1987A is interesting in light of the potential metastability of forming neutron 
stars. On the one hand, on February 23 of 1987 neutrinos were observed from the explosion 
of supernova SN1987A, indicating that a neutron star, not a black hole, was initially present. 
On the other hand, the ever-decreasing optical luminosity of the remnant of SN1987A 
suggests two arguments [64, against the presence of a neutron star. 



First, accretion at the Eddington limit with the usual Thomson electron scattering 
opacity onto a neutron star is already ruled out. Chen and Colgate p6| , however, have 
recently suggested that the opacity appropriate for a neutron star atmosphere has been 
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underestimated by several orders of magnitude. Using the opacity of iron at X-ray photon 
energies, they conclude that the appropriate Eddington limit cannot yet rule out accretion 
onto a neutron star. 

Second, a Crab-like pulsar cannot exist in SN1987A since the emitted magnetic dipole 
radiation would be too large. The magnetic field and/or the spin rate of a neutron star 
remnant must be much less than in the case of the Crab and, therefore, much less than is 
inferred for other young neutron stars. It is possible, however, that although the spin rate 
of a newly formed neutron star is expected to be high, the timescale for the generation of 
a significant magnetic field is greater than 10 years. Unfortunately, this timescale is not 
known with certainty 



The experimental measurements [68] of neutrinos from SN1987A indicated the follow- 
ing: 

• A total binding energy of ~ (0.1 — 0.2)Mq was released, indicating, from Fig. 28, a 
remnant gravitational mass of (1.14 — 1.55)M0. In addition, about half or more of the 
binding energy appears to have been released during the first 2 seconds, in agreement 
with the analysis of the previous section and Fig. 29. The binding energy arguments 
therefore do not discriminate among the various scenarios we have discussed. 

• The average neutrino energy was ~ 10 MeV; to lowest order, this is fixed by the mean 
free path \ V {E V ) in the outer regions of the protoneutron star, and also does not shed 
much light on the internal stellar composition. 

• In spite of the fact that most of the binding energy is released during the initial 
accretion and collapse stage in the first 2 or so seconds after bounce, the neutrino 
signal continued for a period of at least 12 s. This latter timescale may be significant, 
since it is the also about the time required for the neutrinos initially trapped in 
the star to leave. However, counting statistics prevented measurement of a longer 
duration, and this unfortunate coincidence prevents one from distinguishing a model 
in which negatively-charged hadronic matter appears and a black hole forms from 
a less exotic model, in which a neutron star still exists. As we have pointed out, 
the maximum stable mass drops by ~ O.2M when the trapped neutrinos depart if 
negatively charged hadrons are present, be they hyperons, kaons or quarks, which 
could be enough to lead to continued collapse to a black hole. 

Observed neutron stars lie in a very small range of gravitational masses. The smallest 
range that is consistent with all the data |3£| runs from 1.34Mq to 1.44Mq, the latter value 
being the accurate measurement of PSR1913+16. Thielemann, et.al. j7(| and Bethe and 



Brown [65] have estimated the gravitational mass of the compact core of SN1987A to be in 



the range ~ (1.40 — 1.56)M , using arguments based on the observed amounts of ejected 
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56 Ni and/or the total explosion energy. This range extends above the largest accurately 
known value for a neutron star mass, 1.44 M Q , so the possibility exists that the neutron 
star initially produced in SN1987A could be unstable in the cold, deleptonized state. A 
possible scenario is that such a mass could be stabilized initially when the neutrinos were 
trapped, but could become unstable when the neutrino concentration dropped to very small 
values (see Fig. 16, for example). In this case, therefore, SN1987A would have become a 
black hole once it had deleptonized, and no further signal would be expected. 

Note that this scenario for black hole formation in SN1987A is different from that 
originally proposed by Brown, Bruenn and Wheeler f64|] , who suggested that long-term 
accretion on the remnant eventually resulted in the production of a black hole. Although 
recent 2-dimensional hydrodynamical calculations of supernovae 61] suggest that sig- 
nificant accretion ceases when the supernova shock lifts off, no more than 1 second after 
bounce, it remains to be seen whether the reverse shocks generated when the shock reaches 
the low-density hydrogen layers produces significant fallback [71]. 



6.4 Future detections 

A fundamental question is whether or not future neutrino detectors will be able to dis- 
criminate among EOS models. The neutrino signal observed in a terrestrial detector is the 
folding of the emitted neutrino spectrum with the detector characteristics. The latter is 
a combination of the appropriate microphysical cross sections for neutrino scattering and 
absorption together with the efficiency function and fiducial (effective) volume of the de- 
tector. Another important parameter of the detector is its low-energy threshold, or cutoff. 



Burrows, Klein and Gandhi [72] list some properties of present, under construction, and 
future detectors, together with a rough estimate of the numbers of neutrinos that would be 
observed from a supernova located within our Galaxy (assumed to be at a distance of 10 
kpc). Table 17 lists some of the characteristics of the various neutrino telescopes. 

In an optimistic scenario, about 10,000 neutrinos will be seen from a typical galactic 
supernova in a single detector. A crucial question, which has not yet received much atten- 
tion, is whether the statistical uncertainty in a time-dependent signal can be small enough 
to adequately differentiate models. 



Among the interesting features that could be sought are: 

1. Possible cessation of a neutrino signal due to black hole formation. 

2. Possible burst or light curve feature associated with the onset of negatively-charged 
hadrons near the end of deleptonization, whether or not a black hole is formed. 
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3. Identification of the deleptonization/cooling epochs by changes in luminosity evolution 
or in neutrino flavor distribution. 



4. Determination of a radius-mean-free-path correlation from the luminosity decay time 
or the onset of neutrino transparency. 

5. Determination of the neutron star mass from the universal-binding energy-mass re- 
lation. 

To realize the above goals, more information about the characteristics of neutrino 
telescopes must be made widely available. This is especially important in deciphering the 
time evolution of the neutrino signal (see, for example, Lattimer and Yahil f52[| ) even if a 
large number (10,000 or more) of neutrinos are detected in total. 
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Table 1 

Potential model parameters for nuclear matter 



EOS 


K 


A 


B 


B' 


a 


Ci 


c 2 


BPAL1 


120 


75.94 


-30.88 





0.498 


-83.84 


23 


BPAL2 


180 


440.94 


-213.41 





0.927 


-83.84 


23 


BPAL3 


240 


-46.65 


39.45 


0.3 


1.663 


-83.84 


23 


SL1 


120 


3. 706 


-31.155 





0.453 


-41.28 


23 


SL2 


180 


159.47 


-109.04 





0.844 


-41.28 


23 


SL3 


240 


-204.01 


72.704 


0.3 


1.235 


-41.28 


23 



Parameters in Eq. ( |10|) determined by fitting the equilibrium properties of sym- 



metric nuclear matter for some input values of the compression modulus Kq [17] 



All quantities are in MeV, except for the dimensionless a. The finite-range pa- 
rameters Ai = 1.5pp and A2 = 3pp\ 
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Table 2 

Potential model parameters for neutron-rich matter 



EUb 


TV" 


Xq 


£3 


Zi 




BPAL11 




-0.689 


0.577 


-14.00 


16.69 


BPAL12 


120 


-1.361 


-0.244 


-13.91 


16.69 


BPAL13 




-1.903 


-1.056 


-1.83 


5.09 


BPAL21 




0.086 


0.561 


-18.40 


46.27 


BPAL22 


180 


-0.410 


-0.105 


-9.38 


24.05 


BPAL23 




-1.256 


-1.358 


-11.67 


-10.90 


BPAL31 




0.376 


0.246 


-12.23 


-2.98 


BPAL32 


240 


0.927 


-0.227 


-11.51 


8.38 


BPAL33 




1.654 


-1.112 


3.81 


13.16 


SL12 


120 


-3.548 


-0.5 


-13.355 


2.789 


SL22 


180 


-0.410 


-0.105 


9.38 


-4.421 


SL32 


240 


-0.442 


-0.5 


-13.387 


2.917 


in Eq. flll 


); x Q 


and £3 are 


dimensionless, the 


remainin 



are in MeV. For each compression modulus Kq, the three different choices of the 
constants yield the potential part of the symmetry energy that varies approx- 
imately as y/u ,u and 2it 2 /(l + u), respectively, as in the parameterization of 
Ref. [17]. In all cases, the symmetry energy at the nuclear matter equilibrium 



density is taken to be 30 MeV. The notations BPALniri2 and SLnin2 are used 
used to denote different EOSs; n\ refers to different values of Kq, and n<i = 1,2 
and 3 indicate, respectively, a \fu ,u, and 2n 2 /(l + n) dependence of the nuclear 
symmetry potential energy on the density. 
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Table 3 

Pure neutron star properties at finite entropy in the potential models. 



EOS S R ^ P c T c A-10 2 / 

M n 









(km) 




MeV mi" 3 


MeV 




Mq km 2 







1.896 


10.509 


7.344 


546.6 


0.0 




87.60 


BPAL 22 


l 


1.941 


10.987 


6.812 


493.0 


72.9 


2.61 


95.73 




2 


2.093 


12.331 


5.392 


346.4 


138.4 




126.60 







2.020 


10.55 


7.05 


645.2 


0.0 




98.97 


SL22 


1 


2.109 


11.12 


6.37 


565.6 


117.0 


4.31 


113.53 




2 


2.369 


12.59 


4.90 


429.2 


208.2 




163.24 




T c , 


and I 


refer to the radius, central density, 


pressure, 


tempera- 



ture, and moment of inertia of the maximum mass star. The coefficient A (see 
Eq. (p9])) shows the increase in the maximum mass due to thermal effects. 
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Table 4 

Star properties for matter in beta equilibrium at finite entropy 
in the BPAL potential model. 



EOS 


S 


71 /f 

«max 


R 


n c 


Pc 


T c 


A-10 2 








Mq 




n 
















km 




MeV fm~ 3 


MeV 




M Q km z 







1.393 


8.219 


12.656 


843.1 


0.0 




36.49 


BPAL 11 


1 


1.415 


8.540 


11.875 


786.6 


55.9 


1.73 


39.11 




2 


1.489 


9.571 


9.687 


552.4 


97.6 




48.84 







1.454 


8.943 


10.938 


659.5 


0.0 




43.49 


BPAL 12 


1 


1.475 


9.250 


10.156 


621.5 


43.5 


1.48 


46.35 




2 


1.540 


10.219 


8.594 


451.6 


82.4 




56.38 







1.473 


9.446 


10.156 


566.6 


0.0 




47.16 


BPAL 13 


1 


1.495 


9.872 


9.375 


449.7 


39.3 


1.60 


51.00 




2 


1.567 


10.840 


7.812 


376.4 


74.1 




62.67 







1.672 


9.172 


9.687 


766.2 


0.0 




58.83 


BPAL 21 


1 


1.689 


9.437 


9.219 


712.4 


50.5 


1.00 


61.73 




2 


1.739 


10.239 


7.969 


550.9 


88.4 




71.07 




n 

u 


1 799 


Q 791 


8 7^0 


uoo.o 


n n 




uu.ou 


BPAL 22 


1 


1.735 


9.943 


8.437 


564.2 


40.5 


0.79 


68.78 




2 


1.776 


10.630 


7.500 


500.7 


77.1 




77.64 







1.737 


10.104 


8.281 


566.5 


0.0 




70.13 


BPAL 23 


1 


1.752 


10.350 


7.969 


536.4 


36.1 


0.88 


73.26 




2 


1.798 


11.120 


6.875 


409.9 


69.5 




83.63 







1.905 


10.107 


7.734 


652.1 


0.0 




84.79 


BPAL 31 


1 


1.917 


10.280 


7.500 


629.0 


41.3 


0.62 


87.31 




2 


1.952 


10.878 


6.878 


543.5 


78.4 




95.29 







1.933 


10.420 


7.343 


590.2 


0.0 




90.14 


BPAL 32 


1 


1.943 


10.589 


7.138 


577.7 


36.7 


0.53 


92.52 




2 


1.974 


11.136 


6.506 


482.8 


71.5 




100.17 







1.955 


10.797 


7.000 


532.0 


0.0 




95.35 


BPAL 33 


1 


1.966 


11.020 


6.719 


507.0 


33.3 


0.49 


98.57 




2 


1.994 


11.518 


6.198 


454.3 


66.0 




105.72 
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Table 5 

Star properties for matter in beta equilibrium at finite entropy 
in the SL potential model. 



EOS 


S 


M max 

Mq 


R 
km 


n c 
no 


Pc 

MeV fm~ 3 


MeV 


A-10 2 


I 

Mq km 2 







1.7423 


9.145 


9.53 


923.3 


0.0 




62.66 


SL12 


1 


1.7711 


9.526 


8.91 


816.8 


71.4 


1.86 


67.46 




2 


1.8710 


10.61 


7.27 


558.2 


121.0 




85.53 







1.890 


9.840 


8.13 


773.5 


0.0 




80.26 


SL22 


1 


1.920 


10.132 


7.73 


724.9 


64.0 


1.59 


85.15 




2 


2.010 


11.110 


6.56 


500.7 


113.5 




102.60 







2.0971 


10.572 


6.81 


689.9 


0.0 




107.14 


SL32 


1 


2.1211 


10.790 


6.60 


651.6 


56.6 


1.11 


111.60 




2 


2.1901 


11.549 


5.83 


532.2 


103.2 




127.13 
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Table 6 

Star properties for matter with trapped neutrinos (Yl € = 0.4) in beta equilibrium 
at finite entropy in the BPAL potential model. 



EOS 


S 


71 /T 


R 


n c 


Pc 


T c 


A-10 2 








Mq 




n 
















km 




MeV fm~ 3 


MeV 




M Q km z 







1.376 


8.459 


11.875 


766.3 


0.0 




36.48 


BPAL 11 


1 


1.396 


8.705 


11.250 


710.6 


43.3 


1.63 


38.78 




2 


1.465 


9.750 


9.219 


505.3 


78.9 




48.30 







1.394 


8.635 


11.406 


719.0 


0.0 




38.37 


BPAL 12 


1 


1.413 


8.884 


10.781 


663.8 


41.2 


1.59 


40.77 




2 


1.482 


9.919 


8.906 


478.1 


75.7 




50.47 







1.403 


8.782 


11.094 


681.6 


0.0 




39.59 


T >» t- > AT -I l \ 

BPAL 13 


1 


1.422 


9.035 


10.469 


620.5 


39.7 


1.60 


42.11 




2 


1.492 


10.099 


8.594 


447.3 


73.0 




52.37 







1.641 


9.300 


9.219 


705.4 


0.0 




57.38 


T "A T "\ A T r-\ -i 

BPAL 21 


1 


1.656 


9.465 


8.958 


684.7 


38.7 


0.96 


59.40 




2 


1.704 


10.199 


7.917 


549.0 


72.9 




67.80 




n 

u 


± .Uuu 






uoo.u 


n n 






BPAL 22 


1 


1.669 


9.580 


8.750 


664.7 


37.0 


0.93 


61.24 




2 


1.716 


10.330 


7.708 


519.5 


70.3 




69.79 







1.662 


9.505 


8.906 


665.1 


0.0 




60.21 


BPAL 23 


1 


1.676 


9.677 


8.594 


636.5 


35.9 


0.92 


62.45 




2 


1.723 


10.427 


7.604 


509.2 


68.5 




71.06 







1.855 


10.036 


7.656 


644.6 


0.0 




79.65 


BPAL 31 


1 


1.867 


10.198 


7.422 


618.4 


34.1 


0.66 


82.01 




2 


1.904 


10.747 


6.812 


523.9 


66.4 




89.78 







1.862 


10.092 


7.578 


633.2 


0.0 




80.74 


BPAL 32 


1 


1.874 


10.249 


7.344 


610.7 


33.7 


0.64 


83.17 




2 


1.910 


10.820 


6.719 


514.1 


65.3 




90.82 







1.871 


10.203 


7.437 


609.6 


0.0 




82.42 


BPAL 33 


1 


1.883 


10.330 


7.266 


604.6 


32.8 


0.63 


84.44 




2 


1.918 


10.922 


6.625 


503.7 


63.7 




92.35 
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Table 7 

MRHA coupling constants and saturation properties 



M 


M 


K 
(MeV) 


cl 


cl 


r 2 


0.79 


0.66 


354 


180.6 


317.5 


73.5 


1.00 


0.73 


461 


137.7 


215.0 


81.6 


1.25 


0.82 


264 


78.6 


178.6 


90.8 
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Table 8 

Star properties for matter in beta equilibrium at finite entropy 
using relativistic models. 



Mr 


S 




R 

(km) 


n c 


Pc 

MeV fm~ 3 


MeV 


A-10 2 


I 

M Q km 2 


ill 


Mq 


no 







2.532 


12.62 


4.73 


432.2 


0.0 




191.18 


0.79 


1 


2.535 


12.75 


4.63 


416.3 


28.9 


0.10 


193.37 




2 


2.542 


13.00 


4.46 


394.9 


58.3 




197.24 







2.305 


12.02 


5.34 


428.2 


0.0 




151.05 


1.00 


1 


2.311 


12.14 


5.24 


416.6 


27.9 


0.25 


153.19 




2 


2.328 


12.45 


4.98 


387.2 


56.7 




158.82 







1.857 


10.60 


7.29 


484.9 


0.0 




85.63 


1.25 


1 


1.868 


10.72 


7.19 


466.1 


29.1 


0.56 


87.38 




2 


1.899 


11.19 


6.60 


419.6 


58.8 




94.33 







2.005 


10.92 


7.14 


545.8 


0.0 




100.6 


GM 


1 


2.014 


11.08 


6.95 


521.9 


31.6 


0.47 


103.1 




2 


2.044 


11.56 


6.43 


458.2 


62.6 




110.6 



The symbol GM refers to the EOS of Ref. The other three sets are for the 
MRHA model. 
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Table 9 

Star properties for matter with trapped neutrinos (Yj, e = 0.4) in beta equilibrium 
at finite entropy using relativistic models. 



Mr 


S 


M max 


R 


n c 


Pc 




A-10 2 




1V1 




1\/T 

Mq 




no 
















(km) 




MeV fm~ 3 


MeV 




M Q km 2 







2.447 


12.28 


4.86 


460.1 


0.0 




172.77 


0.79 


1 


2.452 


12.29 


4.78 


448.9 


26.4 


0.25 


174.73 




2 


2.471 


12.70 


4.55 


415.7 


52.6 




181.38 







2.222 


11.68 


5.47 


451.2 


0.0 




135.90 


1.00 


1 


2.228 


11.72 


5.42 


448.5 


25.9 


0.40 


137.27 




2 


2.257 


12.16 


5.09 


406.4 


51.5 




145.39 







1.784 


10.28 


7.50 


514.1 


0.0 




76.29 


1.25 


1 


1.791 


10.23 


7.42 


511.4 


27.6 


0.75 


77.16 




2 


1.836 


10.89 


6.81 


448.3 


54.6 




85.25 







1.935 


10.53 


7.41 


595.8 


0.0 




90.13 


GM 


1 


1.946 


10.69 


7.22 


568.3 


30.2 


0.58 


92.48 




2 


1.980 


11.19 


6.67 


496.6 


59.0 




100.1 
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Table 10 

MRHA coupling constants and saturation properties with hyperons. 



M 


M* Nsat 
M 


K 
(MeV) 


r 2 


cl 






0.79 


0.76 


177 


118.7 


258.1 


84.8 


0.660 


1.00 


0.73 


455 


133.1 


210.3 


82.4 


0.658 


1.25 


0.84 


228 


64.9 


174.0 


92.6 


0.679 



In all cases, the ratio of hyperon to nucleon sigma and rho couplings are taken 
to be equal: x a = g aH /9an = x p = g p H/9pn = 0.6. 
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Table 11 

Star properties for matter, including hyperons, in beta equilibrium 
at finite entropy using relativistic models. 



Mr 


S 


M max 


R 

(km) 


n c 


Pc 

MeV fm~ 3 


MeV 


A-10 2 


I 

M Q km 2 


1\/T 
1V1 


1\/T 

Mq 


no 







1.638 


10.62 


7.81 


376.7 


0.0 




67.8 


0.79 


1 


1.639 


10.73 


7.66 


365.6 


19.7 


0.08 


68.6 




2 


1.643 


10.99 


7.39 


353.5 


42.2 




70.0 







1.886 


12.12 


5.64 


248.1 


0.0 




107.4 


1.00 


1 


1.884 


12.16 


5.64 


250.1 


16.9 


-0.11 


106.9 




2 


1.878 


12.32 


5.54 


247.9 


35.8 




106.6 







1.407 


10.43 


8.44 


310.3 


0.0 




51.0 


1.25 


1 


1.414 


10.59 


8.13 


293.1 


18.4 


0.36 


52.5 




2 


1.428 


10.86 


7.81 


283.9 


39.2 




54.5 







1.544 


10.78 


7.66 


311.4 


0.0 




63.2 


GM 


1 


1.551 


10.95 


7.34 


290.5 


19.7 


0.47 


65.2 




2 


1.573 


11.32 


6.88 


269.5 


41.4 




69.2 
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Table 12 

ar properties for matter, including hyperons and trapped neutrinos (Yt e = 0. 
in beta equilibrium at finite entropy using relativistic models. 



Mr 


S 


M max 


R 


n c 


Pc 




A-10 2 


I 


1\/T 
1V1 




M Q 




no 
















(km) 




MeV fin" 3 


MeV 




M km 2 







1.843 


11.03 


6.60 


368.5 


0.0 




89.5 


0.79 


1 


1.837 


11.15 


6.46 


355.4 


17.0 


-0.15 


89.8 




2 


1.831 


11.49 


6.25 


337.7 


36.3 




91.2 







2.066 


12.15 


5.27 


290.2 


0.0 




126.5 


1.00 


1 


2.063 


12.27 


5.19 


281.5 


15.4 


-0.11 


127.1 




2 


2.057 


12.56 


5.03 


266.9 


32.2 




128.8 







1.580 


10.48 


7.66 


355.9 


0.0 




63.1 


1.25 


1 


1.585 


10.61 


7.50 


344.6 


18.0 


0.30 


64.3 




2 


1.599 


11.11 


6.88 


299.7 


36.9 




64.9 







1.768 


11.11 


6.63 


334.8 


0.0 




83.6 


GM 


1 


1.772 


11.21 


6.56 


332.0 


17.5 


0.10 


84.5 




2 


1.776 


11.66 


6.15 


296.7 


37.0 




88.5 
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Table 13 

Critical density ratio, u cr n = n cr i t /no, for kaon condensation in the relativistic 
mean field model, GM, for the neutrino-free and trapped neutrino cases (Yl € = 0.4), 

with and without hyperons. 









Without hyperons 


With hyperons 


a 3 m s (MeV) 


model 


S 


Y v = 


Y Le = 0.4 


Y u = 


Y Le = 0.4 


-134 


chiral 





4.15 


6.38 


9.46 


** 




mes. exch. 





4.54 


7.29 


** 


** 


-222 


chiral 





3.15 


4.35 


4.22 


? 




mes. exch. 





3.59 


5.46 


** 


** 


-310 


chiral 





2.49 


3.15 


2.73 


? 




mes. exch. 





2.86 


4.03 


3.76 


? 



The symbol ** indicates that no condensation takes place for this set of cou- 
plings. 
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Table 14 

Properties of a star, both without and with trapped neutrinos (Yi e = 0.4), which contains neutrons, protons, 
and kaon condensates in beta equilibrium in the relativistic mean field GM model. 











Y„ 


= 






Y Le 


= 0.4 




Model 


a 3 m s 
MeV 


S 


-^max 


R 

(km) 


n c 


I 

M Q km 2 




R 

(km) 


n c 


I 

Mq km 2 


Mq 


no 


Mq 


n 




-134 





1.911 


11.39 


6.38 


98.8 


1.934 


10.56 


7.05 


90.6 


chiral 


-222 





1.781 


9.89 


8.62 


69.4 


1.902 


10.60 


7.05 


88.1 




-310 





1.779 


9.02 


9.78 


63.4 


1.838 


9.94 


8.05 


75.2 


mes. 


-134 





1.950 


11.36 


6.38 


102.2 


1.935 


10.53 


7.05 


90.3 


exch. 




1 


1.971 


11.50 


6.25 


105.1 


1.945 


10.74 


6.84 


92.9 






2 


2.015 


11.79 


5.99 


111.3 


1.977 


11.22 


6.37 


100.1 


mes. 


-222 





1.832 


10.65 


7.50 


81.6 


1.928 


10.65 


6.88 


91.4 


exch. 




1 


1.866 


11.22 


6.72 


90.9 


1.939 


10.78 


6.84 


92.9 






2 


1.946 


12.01 


5.80 


107.8 


1.970 


11.23 


6.38 


99.8 
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Table 15 



Ratios of hyperon-meson to nucleon- meson coupling constants, xw = gm 
where i = a, w or p, and H is a hyperon species. 



Case 




XujA 


x pA 






XpZ 




Xu)=. 


XpE 


1 


0.60 


0.65 


0.60 


0.54 


0.67 


0.67 


0.60 


0.65 


0.60 


2 


0.60 


0.65 


0.60 


0.77 


1.00 


0.67 


0.60 


0.65 


0.60 


3 


0.60 


0.65 


0.60 


0.77 


1.00 


0.67 


0.77 


1.00 


0.67 
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Table 16 



Maximum masses of stars, M max /MQ, with baryonic matter that undergoes a 
phase transition to quark matter without (Y u = 0) and with (Y^ e = 0.4) trapped 
neutrinos. Results are for a mean field model of baryons and a bag model of 
quarks. B denotes the bag pressure in the quark EOS. 





Without hyperons 


With hyperons 


B 

(MeV fm" 3 ) 


Y v = 


Y Le = 0.4 


Y u = 


Y Le = 0.4 


136.6 
150 
200 
250 


1.440 
1.444 
1.493 
1.562 


1.610 
1.616 
1.632 
1.640 


1.434 
1.436 
1.471 
1.506 


1.595 
1.597 
1.597 
1.597 


No quarks 


1.711 


1.645 


1.516 


1.597 
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Table 17 



SUPERNOVA NEUTRINO TELESCOPE CHARACTERISTICS 



Detector 


Total mass (tonnes) 


Composition 


Threshold 


# Events 




(Fiducial Mass (tonnes)) 




(MeV) 


at 10 kpc 


CERENKOV: 










Kill 


3000 


H 2 


5 


370 




(2140) 








Super Kamiokande 


40,000 


H 2 


5 


5500 




(32,000) 








SNO 


1600/1000 


H 2 0/D 2 


5 


780 


SCINTILLATION: 










LVD 


1800 


Kerosene 


5-7 


375 




(1200) 








MACRO 


1000 


"CH 2 " 


10 


240 


Baksan 


330 


"White Spirits" 


10 


70 




(200) 


"CH 2 " 






LSND 


200 


"CH 2 " 


5 


70 


Borexino 


300 


(BO) 3 (OCH 3 ) 3 


~ 0.2 


200 


Caltech 


1000 




2.8 


290 


DRIFT CHAMBER: 










ICARUS 


3600 


40 Ar 


5 


120 


RADIOCHEMICAL: 










Homestake 37 CI 


610 


C 2 Cl4 


0.814 


4 


Homestake 12 I 




Nal 


0.664 


25 


Baksan 37 Cl 


3000 


C 2 Cl4 


0.814 


22 


EXTRAGALACTIC: 










SNBO 


100,000 


CaCOg 




10,000 


JULIA 


40,000 


H 2 




10,000 
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8 Figure Captions 



Fig. 1. The main stages of evolution of a neutron star. Numbers in parentheses refer to 
the stages discussed in the text. 

Fig. 2. Results (left panels for BPAL EOS and right panels for SL interactions) for pure 
neutron matter. Top panels show the neutron effective mass ratio from Eq. (|25| ) and Eq. (|26| ) 
versus the density ratio u = n/riQ, where uq = 0.16 fm -3 is the equilibrium nuclear density. 
The center panels show isentropic pressures, and the bottom panels show star masses versus 
central density ratio at fixed entropy per baryon. 

Fig. 3. Results (left panels for BPAL EOS and right panels for SL interactions) for matter 
in beta equilibrium among n, p, e~ , and \x~ , at an entropy per baryon S = 1. Top panels: 
Individual concentrations Y% = m/rib, where i = n, p, e~ and fjT. Center panels: The 
electron chemical potential \i e = fj,u = fi n — fj, p . Bottom panels: Individual contributions 
to the entropy per baryon. 

Fig. 4. Results (left panels for BPAL EOS and right panels for SL interactions) for matter 
in beta equilibrium among n, p, e~ and fi~ . (See caption to Fig. 2 for further details. 
Proton effective masses are also shown here.) 

Fig. 5. The moment of inertia, /, as a function of density (left panel) and baryonic mass, 
Mb (right panel). The BPAL22 equation of state is employed for fixed values of the entropy 
per baryon. The full dots on the curves indicate the maximum gravitational mass. 

Fig. 6. Results for the BPAL model with trapped neutrinos at entropy per baryon S = 1. 
The upper panel shows individual concentrations, the center panel gives the leptonic chem- 
ical potentials, with fi = fi e — fj, Ue , and the lower panel separates the nucleon and lepton 
contributions to the entropy per baryon. 

Fig. 7. Stellar temperature, T, as a function of the density ratio u for the MRHA model 
with fi r /M=1.25. The full curves are for neutrino free matter, and the dotted curves refer 
to matter with trapped neutrinos. In the upper (lower) panel, the baryons are nucleons 
without (with) hyperons. 

Fig. 8. Results for matter in beta equilibrium among n, p, e~, and fjT in the MRHA model, 
with fi r /M=1.25, at an entropy per baryon 5 = 1. (See caption to Fig. 3 for further details.) 

Fig. 9. Top panel: Nucleon Landau effective mass ratios versus density ratio for the MRHA 
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model (with // r /M=1.25). Middle panel: Isentropic pressures. Bottom panel: Star mass 
versus central density ratio at fixed entropy per baryon. 

Fig. 10. Results for matter in beta equilibrium among n, p, e~, fi~, and trapped neutrinos, 
in the MRHA model, with (i r /M=1.25, at an entropy per baryon S = 1. Shown are the 
individual concentrations (top panel), the leptonic chemical potentials (middle panel), and 
the baryonic and leptonic contributions to the entropy as a function of density. 

Fig. 11. Relative fractions and the electron chemical potential for beta-equilibrated matter 
containing nucleons, hyperons, electrons, and muons in the MRHA model (^ r /M=1.25) at 
zero temperature. 

Fig. 12. Results for matter in beta equilibrium among nucleons, hyperons, electrons, and 
muons in the MRHA model with /i r /M=1.25 at an entropy per baryon 5 = 1. Shown are 
the individual concentrations (top panel), the electron chemical potential (middle panel), 
and the baryonic and leptonic contributions to the entropy as a function of density. 

Fig. 13. Relative fractions and leptonic chemical potentials for beta-equilibrated mat- 
ter containing nucleons, hyperons, electrons, muons, and trapped neutrinos in the MRHA 
model (// r /M=1.25) at zero temperature. Here [i = [i e — Hv e - 

Fig. 14. Results for beta-equilibrated matter containing nucleons, hyperons, electrons, 
muons and trapped neutrinos in the MRHA model (/x r /M=1.25) at an entropy per baryon 
S = 1. Shown are the individual concentrations (top panel), the leptonic chemical po- 
tentials (middle panel), and the baryonic and leptonic contributions to the entropy as a 
function of density. 

Fig. 15. Panel (1): Ratio of gravitational mass Mq to baryonic mass Mb as a function of 
Mb for the GM model. Solid lines are for lepton-rich matter, dashed lines for neutrino-poor 
matter. A dot at the end of a curve indicates matter with hyperons, a star indicates matter 
without hyperons. For the neutrino-poor cases, the entropy per baryon is given next to the 
curves. Panel (2): Gravitational mass as a function of baryonic mass. The symbols are the 
same as in panel (1). 

Fig. 16. Maximum neutron star mass as a function of electron-neutrino fraction Y Ue in the 
GM model for matter with and without hyperons, labeled by npH and np, respectively. 

Fig. 17. Illustrative plot of the kaon energies uj^ in the meson exchange model as a 
function of the density ratio u. Here the baryons are nucleons. The chemical potential fi 



81 



is also shown; the dashed portion of the curve indicates the behavior when kaons are absent. 

Fig. 18. Neutrino-free matter in beta equilibrium among nucleons, (thermal) kaons, elec- 
trons and muons in the GM model, with the meson-exchange formalism as a function of 
temperature, T. Results are shown for three different values of the kaon-nucleon sigma 
term Y, KN . Bottom panel: Critical nucleon density ratio for the onset of kaon condensa- 
tion. Next to bottom panel: The electron chemical potential \i e = [i^ = \x n — fi p = fix- 
Next to top panel: Thermal kaon to baryon ratio at threshold, for kaon condensation. Top 
panel: Proton fraction at threshold. 

Fig. 19. Results for neutrino-free matter in beta equilibrium among nucleons, kaons, elec- 
trons and muons in the GM model for an entropy per baryon S = 1. The top panel shows 
the relative concentrations. The center panel shows the electron chemical potential, and 
the bottom panel shows the contributions to the total entropy from the strongly interacting 
particles and the leptons, respectively. 

Fig. 20. Results for neutrino-trapped matter in beta equilibrium among nucleons, (ther- 
mal) kaons, electrons, and muons in the GM model as a function of temperature, T, with 
three different values of the kaon-nucleon sigma term T, KN . Bottom panel: Critical nu- 
cleon density ratio for the onset of kaon condensation. Next to bottom panel: The electron 
chemical potential fj, e = /j,^ = fj, n — fi p = hk- Next to top panel: Thermal kaon to baryon 
ratio. Top panel: Proton fraction at threshold, for kaon condensation. 

Fig. 21. Results for neutrino-trapped matter in beta equilibrium among nucleons, kaons, 
electrons, and muons in the GM model for an entropy per baryon 5 = 1. The top panel 
shows the relative concentrations. The center panel gives the leptonic chemical potentials, 
with fji = fi e — fj, Ue , and the bottom panel shows the contributions to the total entropy from 
the strongly interacting particles and the leptons, respectively. 

Fig. 22. Results for zero-temperature matter in beta equilibrium among nucleons, hy- 
perons, kaons, electrons, and muons. The chiral model with a^m s = —222 MeV is used 
in conjunction with the mean field GM description of the baryons. (a) Relative fractions 
Yi = n i/(J2B n B)- (b) Baryon Dirac effective masses, the kaon chemical potential fx = fi e , 
and the scalar field a. (c) Baryon scalar densities, (d) Condensate amplitude, 9, in degrees. 

Fig. 23. Particle fractions for model HS81 in conjunction with the kaon meson-exchange 
formalism for different choices of the £ and S coupling constants. Panels (1), (2), and (3) 
correspond to parameter sets 1, 2, and 3 of Table 15, respectively. 
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Fig. 24. Gravitational mass vs. baryonic mass for matter with and without kaons in 
the lepton-rich {Yl s = 0.4) and neutrino-poor stages {Y v = 0). The chiral model with 
azm s = —222 MeV is used in conjunction with the mean field GM description for the 
nucleons. The range in neutron star masses that is metastable during deleptonization is 
indicated. 

Fig. 25. Maximum neutron star mass as a function of electron- neutrino fraction Y Ve for 
matter with and without kaons, labeled by np and npK, respectively. (See caption to Fig. 
21 for further details.) 

Fig. 26. Individual concentrations for matter in beta equilibrium among nucleons, hy- 
perons, quarks, electrons, and muons, employing the mean field GM model in the baryon 
sector and a bag model for the quarks. The top panel shows the neutrino-free case and 
the bottom panel the results with trapped neutrinos. The quark phase cavity pressure 
B = 200 MeV fin" 3 . 

Fig. 27. Quark-hadron phase transition boundaries in beta-equilibrated matter as a func- 
tion of the bag pressure, B. In the top panel, the hadrons are nucleons and, in the lower 
panel, nucleons and hyperons. The onset of a quark-hadron mixed phase occurs at a density 
ratio ui, and a pure quark phase begins at u<i- Also shown is the central density ratio, u c , 
of the maximum mass star. 

Fig. 28. Binding energy versus baryon mass for nucleons-only matter (np), matter with 
nucleons and hyperons (npH), and matter with nucleons and kaons (npK). The stars, dots, 
and triangles mark the maximum mass configurations. The lower envelope is for an EOS 
that is causal above a transition density of n t = 0.3 fm -3 and for the GM EOS below n t . 

Fig. 29. Enclosed gravitational mass versus baryon mass. Two generic compositional 
cases are displayed: nucleons-only matter and matter with hyperons. Solid curves are for 
neutrino-rich matter with Yj, e = 0.4 at an entropy per baryon 5 = 2. Dashed curves refer 
to cold catalyzed neutrino-free matter. 

Fig. 30. Electron concentrations as a function of density for neutrino-free matter with var- 
ious assumptions for the stellar composition, as indicated. The quark phase cavity pressure 
B = 200 MeV fm~ 3 . Arrows indicate central densities of 1.44M stars. Differences of each 
curve from Yl c = 0.4 show the net Y Ue lost at each density during cooling. 
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